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Abstract: The effects of anomalies in high density QCD are striking. We consider a direct 
appHcation of one of these effects, namely topological currents, on the physics of neutron 
stars. All the elements required for topological currents are present in neutron stars: de- 
generate matter, large magnetic fields, and parity violating processes. These conditions 
lead to the creation of vector currents capable of carrying momentum and inducing mag- 
netic fields. We estimate the size of these currents for many representative states of dense 
matter in the neutron star and argue that they could be responsible for the large proper 
motion of neutron stars (kicks), the toroidal magnetic field and finite magnetic helicity 
needed for stability of the poloidal field, and the resolution of the conflict between type-II 
superconductivity and precession. Though these observational effects appear unrelated, 
they likely originate from the same physics — they are all P-odd phenomena that stem from 
a topological current generated by parity violation. 
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1. Introduction 



It is well known that anomalies have important implications for low-energy physics: the 
electromagnetic decay of neutral pions vr*^ — )■ 27 is a textbook example. Anomalies reveal 
intricate relationships between topological objects such as vortices, domain walls, Nambu- 
Goldstone bosons, and gauge fields and often result in very unusual physics. A particularly 
relevant example is the superconducting cosmic string on which an electric current flows 
without dissipation and carries momentum [1]. The effects of anomalies are well established 
and are reviewed in [2]. 

More recently the role of anomalies in QCD has been studied at finite baryon density 
[3, 4] and similar phenomena have been studied in condensed matter systems [5, 6, 7]. Since 
these first steps many other applications of anomalies in dense QCD have been considered: 
an analysis of the axion physics and microscopic derivation of anomalies [8, 9]; studying 
the vortex structure due to the anomalies currents in neutron stars (type-I versus type-II 
superconductivity) [10]; the charge separation effect at the relativistic heavy ion collider 
(RHIC) [11, 12]; magnetism of nuclear and quark matter [13]; anomaly mediated neutrino- 
photon interactions at finite baryon density [14]; the chiral magnetic effect at RHIC [15] 
and many others. 

For this paper we are interested in the result where a topological vector current (along 
with an axial current) is induced in the background of an external magnetic field when 
the chemical potentials of the right and left-handed fermions are not equal, /ir 7^ ^i. As 
argued in [8] this phenomenon depends only on the presence of chiral symmetry and not 
on whether the chiral symmetry is spontaneously broken. Applications of this induced 
vector current have been discussed in neutron star physics [10] and RHIC related physics 
[15]. In [10] it was shown that a non-dissipating vector current running along the magnetic 
flux tubes may change the behaviour of superconductivity from type-II to type-I, even 
though the Landau-Ginzburg parameter k > l/\/2 suggests type-II behaviour. It was also 
argued that the triangular lattice of flux tubes — the Abrikosov lattice — may be completely 
destroyed due to the helical instability in the presence of the induced vector current. This 
would resolve the apparent contradiction with the precession of the neutron stars [16, 17]. 
For another mechanisms that may resolve this contradiction see [18, 19, 20, 21]. 

Our goal is to present a quantitative analysis of the conditions when parity violation 
[^r 7^ ^J'l) occurs in dense stars and a persistent, topological current is induced. We claim 
that topological vector currents do exist in dense stars and we consider applications that 
may help explain many phenomena in dense stars, such as kicks and toroidal fields. 

It is well known that the weak interactions (where parity is strongly violated) play 
a dominant role in neutron star physics. Producing the asymmetry /i^ 7^ for a given 
process is common in the bulk of neutron stars, but we are interested in coherent parity 
violating effects when the asymmetry appears in macroscopically large regions. Sections 
3-7 are devoted to estimating the induced current in different environments that may exist 
in neutron stars, while Section 8 is devoted to possible applications of the topological 
currents. Readers interested in applications may go directly to Section 8, skipping the 
technical aspects of Sections 3-7. 
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If non-dissipating vector currents are induced they can transfer momentum either by 
escaping the star or radiating photons. Our kick engine continues to work even when 
temperature drops well below T <C MeV. This is because it is the chemical potential 
(fii — fir) 7^ that drives the kick, not the temperature. This leads to long sustained kicks, 
rather than short natal kicks. In Section 8 we use a rough calculation to demonstrate that 
our kick mechanism can explain the proper motion of hyperfast pulsars with v > 800 km/s, 
particularly B1508+55 moving at v = 1083]']gQ^ km/s [22]. A detailed calculation can be 
found in [23]. See [24] for a review of kick velocities. 

The kick naturally produces a magnetic field-momentum correlation {P ■ B). Though 
it is believed that the spin and magnetic field are correlated somehow, the sustained nature 
of the kicks allows a rotation- kick correlation (P-^) to form regardless of the angle between 
the star's rotation and magnetic field. A neutron star spins on the order of milliseconds 
and the kick occurs over hundreds of years. This gives rise to a cylindrical symmetry which 
causes the the force of the kick to average along the axis of rotation. The nature of the spin- 
magnetic field correlation can produce kicks that are either aligned or anti-aligned with the 
spin axis. Regardless of its direction, this correlation is a P-odd effect, which indicates it 
must be generated by a P-odd mechanism such as the parity violating processes that power 
our current. In fact, long, sustained kicks are observationally supported by the analysis in 
[25] and the spin-kick correlation is observationally supported in [25, 26, 27]. This can be 
considered indirect support for our proposal. 

Our mechanism is similar to neutrino emission in that we are ejecting a particle from 
the star, but there are some key differences. Neutrino emission is automatically asymmetric 
with respect to the direction of the magnetic field B; however, most neutrino-based kick 
mechanisms have difficulty delivering the produced asymmetry to the surface of the star. 
Only when it can reach the surface can the asymmetry push the star. In most neutrino 
based mechanisms the star must be very hot (T > MeV) for neutrinos to be energetic 
enough to transfer sufficient momentum and the kick occurs in a matter of seconds. But 
at such high temperatures the neutrinos cannot escape the star without interacting and 
washing out the asymmetry. 

We also suggest that a current running along the magnetic flux (the poloidal field) 
may be the source of the toroidal magnetic field and the finite magnetic helicity thought to 
be required in a neutron star [28, 29, 30, 31]. The magnetic helicity, % = j dP'xA ■ B, that 
arises from the linking of toroidal and poloidal magnetic fields is a topological invariant and 
is a P-odd effect that must be generated though parity violation. Our method of gener- 
ating a toroidal field naturally produces magnetic helicity without requiring arguments to 
temporarily break the topological invariance of the helicity. Many attempts at generating 
magnetic helicity rely on instabilities in the magnetic field caused by the star's rotation. 
Such correlations {B ■ Q) are P-even, and though they may generate toroidal fields, they 
cannot be responsible for helicity. 

These applications of topological currents to seemingly unrelated P-odd effects will be 
discussed in greater detail in Section 8 and only comprise a sample of possible applications. 
Remarkably, the effects of topological currents introduced in this paper can be experimen- 
tally tested in terrestrial laboratories, particularly in some condensed matter systems and 
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in RHIC at Brookhaven. 



2. A Big Picture: The Basic Ingredients and Assumptions 

Topological vector currents can couple to charged fermions and become physical, super- 
conducting, electromagnetic currents. Our motivation comes from the possible effects of 
topological currents on the physics of neutron stars. Even though axial currents exist in 
the star they cannot be used as an electromagnetic source and we will not consider them. 
Non-dissipating, induced vector currents at T = have the form [6, 10, 15] 

{j) = {f^l-^^r)^, (2.1) 

where and fii are the chemical potential of the right and left-handed electrons, and $ 
is the magnetic flux. For the current to be nonzero the number of left-handed particles 
must be different than the number of right-handed particles. The weak interaction, which 
strongly violates parity, is a natural source for the required asymmetry and is where we 
focus our attention. 

We will assume that electrons are the only reasonable charge carrier because all other 
charged particles in the star are too heavy. The particles in a neutron star attain equi- 
librium through the weak interaction, which creates predominantly left-handed particles. 
This creates an intrinsic difference in the number of left-handed and right-handed electrons. 
In an infinite system this imbalance would disappear — the average helicity of the electrons 
would be washed out due to the inverse weak P-violating processes. The key is that the 
neutron star is a finite system and electrons are removed from the star by the current 
before they can decay. The asymmetry that created the current is allowed to propagate to 
the surface and not get washed out.^ 

This topological current corresponds to the lowest energy state in the thermodynamic 
equilibrium when m ^ fir is held fixed. In reality there is a tendency for fj,i and fir 
to equilibrate though weak interactions; however, due to the finite size of the system a 
complete equilibrium cannot be achieved. This is analogous to how neutrinos in cold 
stars can leave the system without further interactions, but unlike with neutrinos the 
electron chemical potential does not drop to zero. The induced current only remains non- 
dissipating when the system is degenerate, n ^ T. In the star's crust this condition 
becomes invalid, the current will become dissipative, and the trapped electrons will return 
into the system. A detailed discussion of this region when ~ T is beyond the scope of 
the present work, and will be discussed somewhere else. If the electrons manage to escape 
the star the electron chemical potential will slowly decrease until the current may stop 

^Here and in what follows we neglect all QED re-scattering effects, which are much stronger than weak 
interactions but they are P-even and, therefore, cannot wash out the produced asymmetry. This is discussed 
in detail in Section 3.2. Because of the large magnetic field the electron only travels in the direction of the 
magnetic field while the motion in transverse directions confined to Landau levels. The term "mean free 
path" in this paper implies the weakly interacting P-odd "mean free path'" when a produced asymmetry 
can be washed out. 
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running. Charge neutrality will cause matter to accrete isotropically possibly maintaining 
some of the chemical potential. 

In the following subsections we will discuss the structure and processes of dense matter, 
formally derive the topological current, and discuss in greater detail how the magnitude of 
the current can be estimated. 

2.1 Notation 

We use the convention h = c = ks = 1 unless otherwise stated. We will always denote 
the momentum of particle as pi, the Fermi momentum as ki, and the chemical potential as 
/Xj. When convenience dictates, the subscript i will either be the symbol of the particle or 
a number, which will be labelled on the Feynman diagram. The three momentum will be 
bolded Pi, with a magnitude denoted pi, and the four momentum will have a greek index 

2.2 Topological vector currents 

The purpose of this section is to provide a brief introduction to topological currents and 
explicitly derive the form of the vector current. We will only provide a sketch here and 
refer the reader to [8, 9, 15] for the finer details. 

Consider the vector current = ■ip^'^ip in the presence of a magnetic field pointing 
z-direction, where is a Dirac spinor. Ultimately we are interested in the left and right- 
handed modes and write the current in the Weyl representation, 

j3 = ^pja% - il^y^r . (2.2) 

The Dirac equation in a magnetic field with potential A can be written as 

{H^ + m^)iPr = E'^iPr, (2.3) 

where ipi = ^{E — H)ipr- Eigenvalues of H acting on ipr are labelled e, hence the Dirac 
equation has two solutions, E± = ib\/ + m^. The Dirac spinor can be written entirely in 
terms of the right-handed spinor and its eigenvalues, 

" (v'r)^ " [4(m2 + e2)]V4 

We further break up the operator H into transverse and longitudinal components H = 
p^a^ + H±, where p^ are momentum eigenstates along the magnetic field and H± = {—idi + 
eAijcji is a transverse operator with eigenstates |A). We take the longitudinal direction to 
be periodic in L and take the limit L — )■ oo later. One can prove that the zero modes of 
H_\_ are simultaneously eigenstates of H with eigenvalue e = pscr^. 

The expectation value for the current is found in the usual manner by summing the 
current over all states weighted by the probability of each state. For fermions the probabil- 
ity is given by the Fermi-Dirac distribution. The derivation for the axial current assumes 
that the densities for the left and right-handed modes are equal. Here, in the derivation 



^f!;r^:J»>. (2.4) 
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of the vector current, we consider the possibihty that the densities are different and as- 
sign the left and right-handed modes each their own Dirac distribution, ni and n^, in the 
expectation value, 

E 

= [nr{E+) + nr{E_) - ni{E+) - ni{E_)] V^t^^r , (2.5) 

where we used the spinor definition above and the fact that summing over terms odd in e 
vanish. If we write everything in terms of the eigenstates of H±^ it can been shown that 
only the zero modes of A survive, 

P3 A=0 

Following the standard arguments, the factor (A|(T'^|A) counts the difference in trans- 
verse zero modes travelling parallel to the magnetic field with positive or negative eigenval- 
ues, and A^_. Taking L — )• oo and integrating each dirac distribution gives the number 
density of 1-dimensional left and right-handed fermions, n{T,fi)r — n{T,fi)i. With this all 
taken into account the current can be written as, 

= [niiT,f,) - n,(r,/i)](iV+ - iV_). (2.7) 

The difference in the positive and negative modes travelling along the magnetic field is 
given in physical terms by the index theorem [8], 

^+-^- = ^' (2.8) 
where $ is the magnetic flux. The current is then, 

{j^) = {ni-nr)^, (2.9) 

where nr{T, /i) and ni{T, /x) are one dimensional number densities of left and right-handed 
Dirac fermions. Furthermore, the one-dimensional number density in the massless limit can 
be written in terms of the chemical potential n(T, fi) = fx/n in which case the expression for 
the current is reduced to (2.1). In this case it is purely topological. At T = the density 
one dimension density of massive particles has simple expression n(0,/.i) = \/ f^"^ — m^/vr. 
The mass slightly diminishes the magnitude of the current. 

These currents are not exotic, but are simply statements of the motion of left and 
right-handed particles given their spin alignment in a magnetic field. Formulae (2.5) and 
(2.9) have a very simple physical meaning: to compute the current one should simply 
count the difference between left-handed and right-handed modes in the background of a 
magnetic field, see Figure 1. We assume that the modes of the current couple to electrons. 
In this formulation we only consider the lowest landau level, which means the spin of an 
electron will antialign with the field. Then if there is an excess of left-handed electrons they 
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Figure 1: A left-handed electron placed inside the magnetic flux in a background of electrons with 
no average helicity will be pushed out. The current wants to balance the number of left-handed 
and right-handed particles and acts as a helicity pump. 



must move in the direction of the magnetic field and if there is an excess of right-handed 
electrons they must move against the field. The topological vector current acts as a pump 
to remove the average helicity and it stops pumping once the average helicity is zero again. 

In general the difference in left and right-handed modes is a complicated function of 
many parameters: magnetic field B, chemical potential /i, temperature T, and mass of the 
particle m; but in the chiral limit (m = 0) the expression for the current takes the simple 
form (2.1), has pure topological character, and can be derived from an anomalous effective 
Lagrangian without referring to the dynamics. The current is expressed in terms of one- 
dimensional fermi distributions (2.7) and the physics of two other dimensions is determined 
by Landau levels (the lowest one for m = T = 0). When T ^ and m ^ the current will 
still be induced, but it will not have a simple topological form (2.1). The relevant formula 
in this case becomes much more complicated and is determined by the ratios of a number 
of dimensional parameters mentioned above, see [15] where some limiting cases have been 
studied. 

The current is insensitive to the structure of the magnetic field and will be induced 
if the magnetic field is confined to magnetic fiux tubes or uniformly distributed. The 
current is simply confined to the regions where the magnetic field is present. For our 
purposes it is not essential whether the magnetic field is represented by magnetic fiux 
tubes, as found in type-II superconductors, or by magnetic domains, the typical structure 
for the intermediate state^ argued for in argued for in [10]. The current is strongest in the 
degenerate regions with fi ^ T where the background magnetic field is large and when 
the system becomes a less degenerate (or not degenerate at all) one should expect strong 

^The intermediate state is characterized by alternating domains of superconducting and normal matter 
where the superconducting domains exhibit the Meissner effect, while the normal domains carry the required 
magnetic flux. 
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suppression [8, 15]. For numerical estimates of the effect it is convenient to count number 
of superconducting flux tubes in the entire star and compute the current per unit quantum 
flux $0 = ^ = 2- 10~^ G cm^, where g = 2e is the charge of the proton Cooper pair. 

Finally, we note that two different applications of the induced vector current have 
been previously discussed: a) neutron stars [10] and b) RHIC physics [15]. Our basic 
objective remains the same as in our previous paper [10]; however the main focus in this 
paper will be estimating the magnitude of the induced current itself while our previous 
paper [10] was mainly devoted to the specific application on study of the type-I vs type-II 
superconductivity as a result of induced vector current. 

2.3 Basic w^eak processes 

The primary composition of nuclear matter in a neutron star is constantly being debated. 
Fundamentally a neutron star is made of neutrons with small, equal fractions of protons 
and electrons. In more exotic models hyperons may appear along with pion and kaon 
condensates. In an effort to simplify the discussion we will constrain ourselves to four 
fundamental interactions that describe the majority of cases in dense stars. 

These interactions have been discussed before in the context of cooling a neutron star 
through neutrino emission. The first two are quite closely related — the direct and modified 
Urea processes, 

e' + P^N + Ve (2.10) 
e- +P + N' ^ N' + N + Ve, (2.11) 

which both obey the beta equilibrium condition /ig + MP = ^^N■ The neutrinos created in 
these processes only interact with matter through the weak force, which is so weak that 
the star is transparent to the neutrinos. They leave the star and do not contribute to the 
equilibrium condition. 

The first of these interactions, the direct Urea Process, should be the dominant process 
in normal nuclear matter but it is heavily suppressed because the the particles are unable 
to conserve momentum while remaining on their Fermi surfaces. It is possible for the direct 
process to conserve momentum if the proton fraction in the star is above 1/9 [32], which 
could occur with the appearance of hyperons [33]. In these cases the direct Urea process 
dominates. The modified Urea process is able to conserve momentum through an external 
nucleon and was used in the first neutrino emission calculations [34] . This is the dominant 
electron producing process in normal nuclear matter, but it is very slow and the presence 
of exotic particles introduces processes that create and destroy electrons quicker. 

As the density of matter increases it is likely that kaon [35] and pion condensates will 
appear. We will restrict the discussion to kaon condensates, which appear at much more 
reasonable densities, 3 no, than pion condensates, 300 no, but the phenomenology of dealing 
with the two condensates is almost identical. Electrons are still created and destroyed by 
the previous interactions, but at a much slower rate than an electron decaying"^ into a kaon 

^Here and in what follows the term "electron decay" means the transforming of an electron into a 
neutrino as a result of interactions with surrounding hadrons. 
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and neutrino in the presence of nucleons, 



e- +N ^ {K)- + N + Ve- (2.12) 

This interaction and its inverse process add another equihbrium condition, /ig = /^i^, on 
top of the previously mentioned beta equilibrium. 

The previous three interactions encompass the creation of electrons in almost all pos- 
sible neutron star interiors. The last interaction we consider is the primary source of 
electrons in quark stars. The direct Urea process for quarks are, 

e'+u^d + Ve (2.13) 
e~ + u ^ s + Ve . (2.14) 

Unlike in normal nuclear matter there is no trouble conserving momentum in quark matter. 
The direct process occurs unsuppressed and there is no need to discuss a modified process 
[36]. 

2.4 Simple models for dense matter 

We will review the features of neutron matter that are required for the rest of the paper. 
Most importantly we will summarize the short reviews found in [36, 37, 38, 39] and state 
the values used in the rest of the paper. The simplest model is the non- interacting gas 
model where the ground state of the neutron star, T = 0, is a mixture of neutrons, protons, 
and electrons that is electrically neutral. The baryon density is on the order of nuclear 
density uq = 0.17 fm~'^, which leads to the nucleons and electrons being highly degenerate. 

The particles achieve equilibrium though the Urea processes (2.10). The neutrinos 
produced in these reactions only react weakly in the star and easily leave it at low temper- 
atures. Because of this neutrinos are often assumed to be non-degenerate and the chemical 
potentials in the star satisfy, 

He + fJ-p = fJ-N ■ (2.15) 

Charge neutrality implies that Ue = np, which implies that the Fermi momenta of the 
electrons and protons are equal, kg = kp. This restriction has two important effects. The 
electrons are relativistic and the protons are non-relativistic implying that the chemical 
potential of the proton is much smaller that of the electron. This further implies that the 
density of neutrons is much higher than the electrons and protons, thus a neutron star. 
Assuming that the density of the neutrons is that of nuclear matter then, 

kN = {STT^npf)^/^ « 340 (n/no)^/^ MeV. (2.16) 

Equating the electron and neutron chemical potentials yields, 

ke « ^ 62 (n/no)^/^ MeV. (2.17) 

2m AT 
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In reality there is a correction to the non-interacting model due to the proton being more 
bound than the electron. It is common in literature to assume a value of 

ke ^ 100 {n/nof^^ MeV, (2.18) 

which is the value we will use through out the paper. 

As the density of the star raises above 3no there is the possibility that condensates 
will appear [35]. As the density of the nuclear matter increases, the density of the electrons 
increases to a point where it becomes advantageous to decay into negatively charged kaons 
though the process (2.12). The system reaches equilibrium through the inverse process. 
These processes add an additional equilibrium condition, 

yLe = HK- (2.19) 

Though negatively charged pions are lighter than kaons, they are unlikely to appear until 
much higher densities due to the strong interaction increasing the effective mass [33]. 

At around the same density as kaons appear, light hyperons and muons may appear 
[33] . The existence of hadrons lowers the neutron density, which lowers the ratio of neutrons 
to protons, making it possible for the direct Urea process to proceed unsuppressed, and 
opens up processes such as K ^ e~ + P + which are also not kinematically suppressed. 
These processes occur at about the same rate as the direct Urea processes, so we will take 
the direct Urea process as a reasonable substitute for them. The rate can then be adjusted 
by an integer factor to compensate for additional processes. The appearance of muons has 
no effect of our calculations as it is too heavy to couple to the current. 

When the density gets high enough it is possible that the quarks deconfine — the 
hadrons break down into their constituent quarks. For this paper we will consider only 
the existence of light quarks in the star, which attain equilibrium through the quark Urea 
processes 2.13. The equilibrium conditions are, 

fiu + fJ-e = IJ^d (2.20) 
+ /^e = Ats , (2.21) 

where, as in earlier cases, the neutrinos are not trapped and are not degenerate. The quark 
matter must also be electrically neutral, 

2 11 

Q/e = -Uu- -Ud- -ns-ne = Q. (2.22) 

The simplest models assume that the quark masses are all zero, thus their Fermi momenta 
are equal to their Fermi energies, and the predicted electron density is zero. This, however, 
is not adequate as leptons do exist in the star. Room for the electrons comes from the 
large mass of the strange quark. Though the up and down quarks are relatively light, 
rriu ~ rrid ~ 5 — 10 MeV, the strange quark mass is actually quite large, rus ~ 100 — 300 
MeV. The strange quark is nonrelativistic meaning there will be fewer of them and electrons 
must be present to conserve charge. For this paper will will follow [36] by assuming that 
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the quarks are massless, k. 



•u 



^ kg ^ kd ^ kq, where 




) 



1/3 



kg = (vr^nfo) 



1/3 



MeV 



(2.23) 



where Ub is the baryon number density. For typical densities in the core of the neutron star 
the Fermi momentum is kg ~ 400 MeV. We can approximate the electron Fermi momentum 
using the fraction of electrons to baryons, = rie/rib, which yields, 



The typical value for the electron fraction is = 0.01. 

We have discussed the nature of topological currents and determined the requirements 
for them to exist. We have also given an overview of many types of dense stars, which are 
the only viable places for these current to appear. Now we are ready to discuss the details 
of how the current appears in the star, and what its magnitude would be. 

2.5 Structure of the magnetic field in a neutron star 

The topological current is confined to regions of the neutron star where there is a mag- 
netic field, but the magnetic flux structure inside a neutron star is non-trivial. Directly 
calculating the current would require careful consideration of the helicity of electrons in 
regions with flux, where all electrons created are left-handed, and regions without flux, 
where the helicity can be washed out, and how electrons diffuse from one region to the 
other. In estimating the current in Section 3 we will simplify this by considering what 
would happen if the flux were uniformly and continuously spread throughout the star. The 
magnitude of the current in the entire star depends only on the amount of flux, not its 
structure. The total current leaving the star would be the same if some complex struc- 
ture were present. Bunching flux lines would simply mean there are smaller regions with 
stronger magnetic fields. However, it is important to be aware of this non-trivial structure 
as it will be considered later in Section 8 when we discuss applications of the current. 

A typical neutron star has a large magnetic field, B ~ 10^'^ G. Once it has cooled, 
the star has a temperature T ~ 10^ K that is cold compared to its Fermi energy fj, ~ 
100 MeV ~ 10^^ K. Because of this the the protons are likely superconducting and the 
neutrons form a superfluid. The magnetic field is large enough that it is favourable for 
the flux to penetrate the superconductor, rather than being completely expelled. The 
Meissner effect forces the flux to bundle into into either type-II vortices or type-I domains, 
which is often called the intermediate state. It is generally believed that the protons 
form a type-II superconductor. The Landau-Ginzburg parameter k = A/^ determines 
the type of superconductivity. Typically in a neutron star the London penetration depth 
is A ~ 120 fm and the coherence length of the proton superconductor is ^ ~ 30 fm. 
This creates a ratio k > which indicates type-II superconductivity. For a type-II 

superconductor the magnetic field will penetrate the star by destroying narrow regions of 
superconductivity that each carry a single quantum of flux, $o = 27r/g = 2 • 10~^ Gcm^. 
But there are problems with this picture [16, 17]. It is possible that the system behaves as 



ke = (SFe) 



/3l 



(2.24) 
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type-I superconductor even though the Landau-Ginzburg parameter would suggest type-II 
behaviour [10]. The argument in [10] relies on the electromagnetic interaction between 
current carrying vortices rather than altering the value of the Landau-Ginzburg parameter 

K. 

If the intermediate state is realized in neutron stars the magnetic field distribution will 
be again non-uniform, but the structure would be quite different. The intermediate state 
is characterized by alternating domains of superconducting and normal matter where the 
superconducting domains exhibit the Meissner effect, while the normal domains carry the 
required magnetic flux. The pattern of these domains is strongly related to the geometry of 
the problem, see [10] for details. While precise calculations are required for understanding 
of the magnetic structure in this case one can give the following estimation for typical size 
of a domain as suggested in [10, 18] 

ar^loVnX, (2.25) 

where i? is a typical external size identified with a neutron star core {R ~ 10 km), while A is 
a typical microscopical scale of the problem. Numerically a ~ 10~^ cm, which implies that 
a typical domain can accommodate about 10^ neutron vortices separated by a distance 10~^ 
cm. While the field distribution for the intermediate state and the type-II superconductor 
are very different one should anticipate that the ratio of normal to superconducting regions 
are the same. 

Regardless of the flux structure, there are two regions of the neutron star to consider — 
those with magnetic flux and those without. The total units of quantum flux can be 
estimated as 

N,^ "^^10^^3,2, Bu^i^r^). (2.26) 



$0 '^"VlOi^G, 

The region that a single unit of flux occupies has a radius equal to the London penetration 
depth of the field A ~ 100 fm. This is multiplied by the number of vortices A^v to get the 
total area. If we take a slice of the neutron star perpendicular to the magnetic field we find 
that the ratio of the area occupied by flux tubes is much smaller than the area occupied 
by the void, 

^vortices _ A^vVtA^ TtX'^B 3 



~ 10~^ • B12 . (2.27) 



This suppression essentially reflects the difference between typical magnetic fleld B ~ 10^^ 
G and the critical magnetic fleld Bd ~ 10^^ G when the superconductivity is destroyed. 



3. Estimating Non-dissipating Currents in Dense Matter 

There are three requirements for topological vector currents to be present: an imbalance in 
the number of left and right-handed particles ni ^ rir, degenerate matter fj,^ T, and the 
presence of the background magnetic fleld B ^ 0. All of these are present in neutron and 
quark stars. The weak interactions, which the star attains equilibrium through, violate 
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parity; particles created in this environment are primarily left-handed, see the Appendix 
for a quantitative estimates. As discussed the interior of the star is dense, /ig = 100 MeV, 
and cold, T = 0.1 MeV, such that the degeneracy condition ^ T is met, and the star is 
known to have a huge magnetic field, B ~ 10^^ G. 

All three criteria are met but there is a subtlety to consider. In an infinite system, 
a system large enough to allow the electron to decay as discussed in section 2.3, any 
asymmetry in left and right-handed electrons created by the weak interaction would be 
washed out; the creation and annihilation rates of the left-handed particles are the same. 
Though many more left-handed electrons are created, they are also destroyed such faster 
than the right-handed and no asymmetry builds. This is similar to the argument found in 
[40] — there is no asymmetry in equilibrium in an infinitely large system. Unlike in [40], we 
are interested in low temperatures where the neutron star is a finite system with respect to 
the weak interaction — see Figure 2. The weak interaction can no longer remove the parity 
introduced into the system. As a result a current forms to remove it. We will first explain 
how the topological currents are induced. We will then discuss why the induced current is 
not washed out by fast quantum electrodynamic (QED) processes, but can only be washed 
out by slow weak interactions. The reason is the P-odd nature of the phenomenology: an 
effect induced by P-odd forces, can only be washed out by P-odd forces. 

mean free path 

in infinite volume ie 




Figure 2: The neutron star is a finite system — the mean free path with respect to the weak 
interaction is much larger than the radius of the neutron star €weak ^ R- Electrons leave the system 
before they decay and contribute towards the current. The current flows because the system is not 
in equilibrium. 



3.1 How topological currents are induced 

If we assume that the magnetic field is uniformly distributed (as discussed in Section 2.5) 
then every electron created by a P-odd process in the star potentially contributes to the 
current. As discussed in Section 2.2 we can calculate the current by counting the number 
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of left-handed electrons minus the number of right-handed electrons created in the star. 
Unlike Section 2.2 there is now a large Fermi momentum that opens up many Landau levels, 
but only the lowest Landau level contributes to the current. The current arises because 
the lowest Landau level has a single spin state, electrons created in that spin state are 
primarily left-handed, and this helicity state propagates out of the star, preserved through 
through countless QED interactions, before the inverse weak process can remove it. These 
helicity states propagating out of the star make up the current. 

We explicitly see how both the magnetic field and parity violation are necessary for 
the current. The spin degeneracy of the lowest Landau level is one, while all other Landau 
levels have spin degeneracy two. This implies that the produced P-asymmetry in the 
polarization (A) = (c? • P /\P\) is not translated into P-asymmetry in the correlation of the 
momentum and magnetic field {B ■ P) for all Landau levels except the lowest one. Even 
though the P-asymmetry is present in the higher levels the spin degeneracy allows particles 
with the same polarization are allowed to travel in opposite directions, which results in zero 
current. The single spin degeneracy of the lowest Landau level means that any longitudinal 
polarization will result in more particles moving one way than the other, thus a current. 
The correlation between spin and the magnetic field is a P-even effect (c? • B) and together 
with the P-odd correlation between spin and momentum {a ■ P) it produces the P-odd 
asymmetry {B ■ P) we are interested in. 

Equation (2.9) tells us how to proceed. The current is found by counting the difference 
in the creation rate of the left and right-handed electrons created in the lowest Landau 
level. This results in a current of the form, 

w 

(j) =i^asym(5,/U,T)- -Fstar, (3.1) 

where is the transition rate per unit volume assuming the magnetic field is uniformly 
distributed, Ps^symiB, fj,,T) accounts for both the polarization of electrons created in a 
magnetic field and the suppression due to Landau levels, and Kitar is the volume of the 
region of degeneracy, which we assume is about the same size as the star itself. 

The extreme degeneracy of the electrons means that electrons are uniformly distributed 
throughout the Landau levels. As discussed in Section 4.2 the number of Landau levels is 
nma.x = (-^e ~ rn^) Bc/2m'^B and the ratio of the lowest Landau level to the total number of 
states is ~ m^B / fi^Bc, which is supported by the analysis in [41]. Therefore, we estimate 
P-odd asymmetry factor as 

T, B) . -(A) . (II) ^ 2 . 10- (I) (^) . (3.2) 

where the polarization (A) is calculated in the appendix and is given by (A. 20). Numerically 
we see that a large chemical potential means that the number of states in the lowest Landau 
level is small compared to the total number of states. We have calculated the average 
polarization not just for the lowest level, but for electron created in all levels. Because we 
are looking at an average current in a large system it is likely that polarization states other 
than those in the lowest landau level are the ones to propagate. What is important is that 
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there are more spin-down states than spin-up states and there is an average polarization. 
The density dependence in (3.2) reflects that found in a neutron star; in the case of quark 
stars we use the appropriate electron chemical potential and density dependence described 
in Section 2.4. The suppression factor is identical to (4.34) but appears because we are 
only considering the Landau levels that contribute to the current. This is different than in 
(4.34) where the factor appears because only one Landau level exists in the system. 

For non-trivial vortex structures it is convenient to determine the current per unit 
fundamental flux by dividing by the total number of flux tubes. In a type-II structure this 
would be the current that runs along a single vortex; in a uniformly distributed environment 
it is simply a convenient normalization, 

(j)=^'asym(B,/^,T).^^. (3.3) 

The current in a type-I domain is found by multiplying (3.3) by the units of flux trapped 
inside the domain. 

The expression for the induced current, (3.1) and (3.3), naively has a different form 
than previously discussed (2.7), but in fact is precisely the same induced current with the 
same physical meaning. The separation of transverse and longitudinal degrees of freedom 
(with respect to magnetic field Bz) that is explicit in (2.7) is hidden now in the formula 
for Pasym(^, A*, T) where the Landau levels (transverse degrees of freedom) are treated sep- 
arately from longitudinal motion. See the Appendix for details. The longitudinal degrees 
of freedom in eq. (2.7) are represented at T = by the one-dimensional number density 
~ ^/vr, which is the correct expression when the problem is treated as a grand-canonical 
ensemble with and /i/j constant due to the infinitely large bath surrounding the system. 
In our case the neutron star is a finite system at T 7^ where particles are continuously 
injected at a rate w. Along with Pasym(-B, /x, T) this describes the resulting asymmetric 
number density as a function of external parameters T,fj,,B. And as they should, both 
equations (2.7) and (3.1) have units of current: number of particles per unit time. 

The topological vector current^ arises specifically because the system is no longer in 
equilibrium with respect to the weak interaction and a small asymmetry (^l — hr) 7^ 
appears. The current is a steady state^, but constantly requires new polarized electrons to 
be created in the magnetic field background, then pushed out of the system, into the crust 
or into space. Because of this steady state the rate w/Q is calculated when the electron 
chemical potential is constant. 

The electron chemical potential does not necessarily go to zero as it does for neutrino 
emission. Only a small fraction of electrons equal to the number in the lowest Landau 
level actually participate in the current. When electrons leave the region of degeneracy, 
where the non-dissipating current is produced, the current loses its quantum properties and 

*The axial current will always be induced even when /ii — fir, but it will not be coupled to the elec- 
tromagnetic field, and cannot carry the momentum. The physical consequences of this axial topological 
current might be quite interesting, but shall not be discussed in the present paper. 

^We assume that the variation of (fiL — Mfl) is adiabatically slow process, with a typical variation time 
to be much longer than any other time scales of the problem. It allows us to treat the system as being in 
the equilibrium when (^l — 1-i.r) is assumed to be a fixed parameter. 
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becomes a normal dissipating current capable of transferring momentum, emitting photons, 
etc. These electrons may get trapped in the crust and diffuse back in, maintaining fig. If 
the electrons manage to escape, the chemical potential //g will slowly decrease over time. As 
electrons escape the star becomes positively charged and stars to accrete matter, which may 
diffuse back into the star. Eventually, when becomes sufficiently small, the production 
of the induced current stops. This is precisely the moment when the neutron star kick 
engine stops. As we mentioned previously it is the chemical potential that fuels the kick, 
not the temperature. The engine stops when chemical equilibrium is achieved, not when 
the temperature drops which happens much earlier. 

3.2 P-odd effects, QED, polarization, and thermodynamics 

An important aspect of the topological current is that it is neither a zero temperature effect 
nor a chiral effect. It persists as long as parity, an asymmetry in left and right-handed 
electrons, is present in the system. Indeed, while the topological currents originally were 
computed at T = for exactly massless fermions [4], it was later shown [8] that the effect 
persists for T ^ 0, rrif. ^ 0. The important lesson learned from these calculations is that the 
effect is not washed out, even when a nonzero temperature is introduced and an arbitrarily 
large number of collisions between electrons occur to maintain thermodynamic equilibrium. 
The crucial point for our mechanism is that P-odd effects are not washed out by fast QED 
processes. 

It is important to remember there are three scales in our problem: the mean free path 
of quantum electrodynamic processes ^qed, the mean free path of weak interactions ^weak, 
and the radius of the neutron star R. For the temperatures we are interested the scales are 
ordered as ^weak ^ ^qed- We see that with respect to QED interactions the neutron 
star is considered to be an infinite system, but with respect to weak interactions it is finite. 

Quantum electrodynamic processes cannot wash out polarization even though massive 
particles introduce processes that can flip helicity, even when these processes occur ther- 
mally. The key here is the unitarity of interactions; processes go forwards and backwards 
at the same rate. Though helicity is lost through these interactions, it is also created. This 
cancellation due to unitarity is aptly illustrated in [40], where the authors consider the 
weak interactions inside a star at high temperatures. Just as in [40] we can use thermal 
equilibrium and unitarity to show how the current propagates. 

Consider electron-proton scattering where neither particle is confined to Landau levels. 
The change in the number of electrons with a given helicity n^'^^(q • cr, t) can be written as. 



where q and p denote the momentum of the electron and proton, cr and s denote the 
electron and proton spin, and W{q' ■ cr' , p' ■ s'|q • <t, p • s) is the probability of scattering 
|p(q, s)e(q, cr)) — )• |p(q', s')e(q', cr')) per unit time per unit phase volume. These are all the 
processes that change the helicity of the electron. 




(3.4) 



[S W{ci • cr, p • s|q • cr', p' • s') - S' VF(q' • cr', p' • s'|q • cr, p • s)] , 
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There is a statistical factor S to account for the Fermi distribution of the initial states 
and the Pauli blocking of the final particles. In the massive limit the chemical potential 
enters the Lagrangian by a term fiN, where the particle number is the charge associated 
with the chemical potential. The mass term literally mixes up the left and right spinors in 
the Lagrangian. Because of this only one value can enter the statistics of the the problem. 
But the chemical potential can be written as fig = + ^^r^ where the left and right-handed 
chemical potentials are really the contribution towards the chemical potential rather than 
chemical potentials in their own right. The current is given by equation (2.9) where we 
calculate the difference in the numbers of left and right-handed electrons. It is common 
to work in the chiral limit where there is no mass term and it is possible in the Weyl 
representation to introduce two separate chemical potentials, one for the left-handed and 
one for the right-handed spinors. There are then two charges that enter the statistics and 
there are two separate Fermi surfaces. The current could then be determined by looking 
at the difference in these chemical potentials. But in the chiral limit there are no helicity 
flipping amplitudes and the imbalance in left and right-handed chemical potentials created 
by the weak interaction would never be washed out. The purpose of this section is to show 
that the current does not get washed out when helicity flipping amplitudes exist so we 
choose the massive limit where the electrons feel only a single chemical potential regardless 
of their helicity. 

With this in mind the statistical factor written explicitly is, 



The factor S' is similar but has the primed and unprimed variables swapped. We have 
chosen degenerate functions to illustrate that the average polarization of the electrons along 
a specific direction (in our case along the magnetic field) is constant, and cannot be washed 
out by QED interactions, see eq. (3.13) below. We have made no assumptions about the 
initial and final chemical potentials and assume that they are different. 
The general statement of thermal equilibrium is 



Thermal equilibrium is attained with respect to QED interactions because the mean free 
path is much smaller than the size of the neutron star. This in not the case with the weak 
interactions, which have a mean free path much larger than that of the star. This non- 
equilibrium of the weak interaction generates the parity asymmetry required for the current 
and is discussed in detail in Section 3.1. Given our specific distributions, the statement of 
thermal equilibrium has two solutions. 



S 



1111 



(3.5) 



s = s' 



(3.6) 



E{p'e) - fJ-'e = E{pe) - fJ-e and E{p'p) - Hp = E{pp) - fip 



(3.7) 



or 



E{pe) - fJ-e = E{pp) - fip and E{p^) - /x^ = E{pp) - Hp . 



(3.8) 
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Substituting either of these solutions into the conservation of energy equation for the 
interaction E{pe) + E{pp) = E{p'^) + E{p'p) yields, 

Pe + f^p = + l^'p , (3.9) 

which is the condition for chemical equilibrium. Particles can only scatter from their Fermi 
surface onto another Fermi surface. 

The key element of the argument is the unitarity of the interaction, 

1= E / ^^VF(q.«T,p.s|q'.^',p'.s') (3.10) 
= Y.I ^^W^(q'-'T',p'-s'|q.^,p.s), (3.11) 

<T',S' % 

This says that every initial state scatters into a final state and that every final state scatters 
into an initial state. 

Using (3.6) the electron and proton distributions can be factored out and using (3.10) 
we find that the forward and reverse rates cancel each other, causing the right hand side of 
(3.4) to vanish. We have found that the number of electrons of a given helicity is constant, 

|n(^)(q.a,t)=0. (3.12) 

This is the result of detailed balance, which ensures that there can be no asymmetry in the 
creation and annihilation rates of particles in thermal equilibrium. Given this spectrum of 
static solutions, the average polarization of the electrons along a specific direction (in our 
case along the magnetic field) is constant, 

/^^q [n(^)(+|q-tT|)-n(^)(-|q-tT|)] 

(A) = -pT:^ — r ^ w , TT , , TTT = constant. (3.13) 

/d3q [n{<=)(+|q • cr|) + n(^)(-|q • cr|)] 

The physical meaning of equation (3.13) is very simple: if we start with no average polar- 
ization, none develops; if we start with an averall polarization, it cannot be washed out. 
This interaction by interaction proof provides the dirty details for a very simple argument: 
a system that starts with some kind of parity violation (i.e., different numbers of left and 
right-handed electrons) cannot have that parity violation removed by QED because QED 
does not violate parity. 

An example of this principle can be found in the experimental set up of a beam of 
electrons. Consider an unpolarized beam of electrons travelling along the z-direction. While 
there are many electromagnetic interactions, the result (3.13) states that the longitudinal 
polarization Vz remains zero as time passes. Why does this happen? Consider the collision 
of two particles in the beam with initial momenta ki and /c2 and zero total polarization 
along the z-axis. After the collision polarization can be induced in the transverse direction 
V ■ ki X /c2, but not in the longitudinal directions, V ■ ki or V ■ k2, as this would contradict 
the fundamental parity symmetry of the QED Lagrangian. In particular, suppose {ki, ^2) 
are in the xz-plane. The polarization can be only induced in the positive y-direction. It is 
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clear that in the thermodynamical equihbrium there will be another pair of particles which 
produce a polarization in the negative y-direction such that total polarization remains zero. 
This simple example explains what equation (3.13) states: if longitudinal polarization of 
the beam was initially zero, it will remain zero in spite of the fact that each given process 
may induce polarization. 

So far the discussion has considered only particle-particle QED scattering and not scat- 
tering off a background field. It is tempting to think that the helicity flipping amplitudes 
of an electron scattering off a semi-classical background field will wash out an asymmetry. 
Such arguments neglect that polarization is given to the background field in the form of 
operators such as the magnetic helicity {A ■ B) and higher order operators that are not 
invariant under P transformations. These operators are the manifestation of the parity 
lost by electrons to the background field, and is the parity responsible for creating the 
toroidal magnetic fields discussed in Section 8. There is a point where the background 
field, which now contains P-odd parity configurations in the form of magnetic helicity, will 
start to give parity back to the electrons through scattering. This steady state is once 
again achieved through unitarity and thermal equilibrium. This system of topological cur- 
rents combined with an E&M field carrying magnetic helicity is complicated and it is not 
our goal to present a complete description of how these P-odd effects transform from one 
form to another. This would require us to analyse a system comprising of time dependent 
Maxwell equations with non-static sources. We explicitly calculate the current when these 
interactions with the external magnetic field are effectively turned off. However, from the 
P-odd invariance of QED we know that P-odd effects will stay even when these interactions 
with the external magnetic field are turned on. 

This has a physical analogue in beam physics as well. As discussed, a beam with 
no longitudinal polarization that has only particle-particle QED interactions will remain 
unpolarized. However, a beam with longitudinal polarization can be created if a magnetic 
field configuration with nonzero helicity is applied, which is the standard technique used 
to produce a longitudinally polarized beam. As this polarized beam propagates it cannot 
lose its polarization due to the internal particle-particle QED interactions. If it did, time 
reversal symmetry, which is unitarity, would tell us that a beam with no longitudinal 
polarization could spontaneously polarize due to particle-particle interactions. This is 
untrue and the polarization stays in the beam. Furthermore, time reversal symmetry 
applied to the longitudinally polarized beam tells us that when the beam encounters a 
magnetic field, the field may reconfigure itself to absorb some P-odd magnetic helicity from 
the beam. This is exactly how the current in a neutron star imparts magnetic helicity to 
the field. To conclude, the entire P-odd configuration which includes longitudinal electron 
polarization and magnetic helicity (and higher order P-odd operators) cannot be washed 
out by QED interactions. 

The final consideration is that some collective bosonic mode such as heat may destroy 
the current. To destroy the current these dissipative modes must carry the helicity out of 
the star faster than the current does. It does not imply that P-odd effects are destroyed. 
Rather it means that P-odd configurations in principle may leave the system. We can 
account for these dissipative modes by assuming they are carried out of the star by photons 



-19- 



or phonons. The helicity modes that the current carries out of the star are subject to direct 
walk and leave the star close to the speed of light while a photon is subject to a random 
walk. The number of steps in a random walk is the square of the number of steps in a 
direct walk, so it takes much longer for a photon to escape than the current to carry helicity 
out. Some tiny fraction of the helicity manages to escape this way, but the current remains 
intact. Thermal cooling of the star is also not the dominant cooling mechanism of the star 
while the current is active, neutrino emission is. The processes that create the neutrinos 
are the exact processes responsible for the current. We can quantify when the star is too 
hot for the current to propagate by finding the temperature where ie < R- This means 
that the current can no longer carrying the helicity out of the system and does not flow. 
Essentially this violates the conditions formulated above for the current to be induced. 

The statement that P-odd effects, if produced, cannot be destroyed or washed out by 
conventional QED processes is correct for any system, including condensed matter systems. 
However, there is a crucial difference between our discussion of the current and conventional 
condensed matter systems, where it is known that the induced spin polarization inevitably 
relaxes even though QED preserves parity. In any condensed matter system the relaxation 
process takes relatively short period of time as the heat can easily leave the system. This 
heat is actually represented by long wavelength photons (microwaves ) that are partly 
polarized and can easily leave the system. The situation is drastically different for neutron 
stars or quark stars (due to the very large density) when it takes very long period of 
time before a photon can reach the surface of the star to have a chance to escape. In 
principle, the P-odd effects in neutron/quark stars will also inevitably relax when no new 
polarization is pumped in. However, the time scale for this to happen is much longer than 
in condensed matter systems. This is analogous to the evolution of the magnetic field in 
neutron stars trapped in a type I superconductor. It is known that the magnetic field 
should be expelled from the bulk of type-I superconductors. However, it takes a very long 
time (much longer than the life time of the universe) before the magnetic field is actually 
expelled from neutron stars. 

Nowhere in the arguments above have we assumed (or implied) that our state is a 
pure quantum state with definite parity; on the contrary, our arguments are based on 
thermodynamics and a density matrix formalism where the polarization for a mixed electron 
state is defined as the sum over all particles V = Qn'Pn where Qn is a probability to find 
the n— th particle with polarization Vn with normalization Yin 9"' ~ ^■ 

For the sake of curiosity there is an example of a strongly interacting, many body 
system where a P-odd configuration could be produced, but nevertheless is not washed out 
by very fast strong interactions, which is very similar to the case we consider. Specifically, 
we have in mind the charge separation effect in heavy ion collisions [11, 12] where P and CP 
odd effect survives (and even have been experimentally observed) in spite of the fact that 
the system is in thermodynamical equilibrium with respect to strong interactions. There 
is a simple argument why a P-odd effect is not washed out by the strong interactions: 
it is an invariance of QCD with respect to P and CP symmetries. The analogy between 
our topological current and the P-odd effect in heavy ion collisions is even deeper than it 
appears: both effects in fact originate from the same anomaly[ll, 12]. 
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4. Nuclear Matter: The Direct Urea Process 



We want to determine the mean free path of the weak interaction of an electron travelUng 
inside a neutron star and the rate at which electrons are created. To do so we will estimate 
the transition rate following the standard techniques from [34]. Estimating the mean 
free path allows us to determine whether the helicity built from the weak interaction is 
washed out or if the asymmetry can escape the star^. Similar calculations have been 
done only for neutrinos as the electron's mean free path is assumed to be much shorter 
due to electromagnetic interactions. However, as we argued above, the electromagnetic 
interactions do not wash out P-odd asymmetry and the electrons are allowed to propagate 
due to the non-dissipating, topological vector current. In order to find the mean free path 
we will consider the direct Urea process, equation (2.10), given by the Feynman diagram 
in Figure 3. 




Figure 3: The direct Urea proeess. 



Given the order of the diagram this should be to dominant process in normal nuclear 
matter, but it is not. It is suppressed because the particles taking part are unable to 
conserve momentum. If all interacting particles lie on their Fermi surface then, 

Pe+pp-pN <-T , (4.1) 

where T is approximately the energy of the neutrino. In order for the process to conserve 
momentum the initial electron and proton, or the final neutron, must be far from their 
Fermi surface in a region with almost no particle occupation. Forcing this inequality to 
become an equality introduces a suppression of order ~ q-^n/t ^ When the proton fraction 
is sufficiently large, Xp > 1/9, this interaction proceeds unfettered [32]. The transition rate 
is given by, 

w = I d^pid^p2d^p3d^pi S (2^)4 n 6{pf - p{) , (4.2) 



Here and in what follows we do not assume that an electron physically escapes a star: we use the term 
"escape" to emphasize that the electron can leave the region of degeneracy without being re-scattered by 
dense surrounding matter. The fate of the moving electrons when they enter the non-degenerate region 
/ie ~ T from deep degenerate region depends on specific properties of matter with ~ T. In this region 
the current becomes dissipating, and the electrons may transfer their energy/momentum to the surrounding 
dense environment. This subject is beyond the interests of the present work, and shall not be discussed 
here. 
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where is the volume of the phase space, the pi are respectively the electron, proton, 
neutron, and neutrino momentum, 5 is a statistical factor that takes into account the 
Fermi blocking, and S the scattering matrix. 

The Fermi blocking factors limit the phase space in which the initial particles can 
exist and the final particles can be created. Particles that exist at the beginning of the 
reaction, such as the electron and proton in this case, are given a factor equal to the Fermi 
distribution, 

= 1 + e-(i?,~/.,)/T ' (^-^^ 

which tells us where particles exist. The particles to be created, such as the neutron, are 
given a factor of one minus the Fermi function, 

^-fi= 1 + e(i?,-M.)/T ' ^^-^^ 

which restricts the phase space that particles can be created in. There is no Fermi blocking 
term associated with the neutrino because they leave the star without interacting and do 
not become dense enough to form a Fermi surface. All of these blocking factors make the 
statistical factor S. 

We use the standard four Fermi scattering matrix element, 



|5|^ = §(l + 3Ci) 



(4.5) 



1 + 3C^ EiEi 

where Ca = 1-26 is the Gamow- Teller coupling and Gp = 1-17 x 10^^^ MeV~^. 

Following [34] we separate the angular and radial integrals such that the transition 
rate becomes. 



where, 



P = jpldpipldp2pldp3pldpiS5{Ef-Ei), (4.7) 
Q = j dnida2da^dni5^^\pf -Pi) . (4.8) 

We start with the Q integral. The contribution from the momentum of the neutrino is 
small compared to the rest so we neglect it in the (5-function. This allows us to take an 
integral over all angles causing the term with angular dependence to vanish. As in [34] the 
angular integrals become, 

Q-^. (4.9) 

We can now do the PQ integral, 

PQ = j pidpip2dp2P3dp3pldpiS6{Ef - Ei) . (4.10) 
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The first step is to change variables from energy to momentum. In doing so we will 
approximate the energy by the value at the Fermi surface as these are the particles most 
likely to participate in the interaction. For the proton and neutron, which are nearly 
non-relativistic, we will approximate the energy by the effective mass, ~ 0.8mN and 
m*p ~ O.Srn-Tv [38]. For the electrons, which are highly relativistic, the chemical potential 
is equal to the Fermi momentum, kg = fig. These factors can be pulled out of the integral. 
Doing the transformation yields, 

pidpi = i^iedEi (4.11) 
P2dpi = m*pdE2 (4.12) 
psdpi = m*^dE3 (4.13) 
pldpi = EldEi. (4.14) 

Performing the neutrino integral, di?4, over the (^-function leaves, 

PQ = ^^^^,m*^m*p j dEi dE2 dE^ {Ei + E2 - E^f S . (4.15) 

The next step is to make the integral dimensionless making the substitutions, 

xi = {El- iie)/T (4.16) 

X2 = {E2- ^ip)/T (4.17) 

X3 = -(^3-/^^)/T, (4.18) 
such that the statistical factor accounting for the Fermi blocking becomes, 



n ("f) 



s 

-«■-«■ l + e 

1=1,2,3 

The Jacobian of these transformations introduces a factor of T for each measure. Also, a 
factor of comes from the {Ei + E2 — E^)"^ term. The chemical potentials introduced 
all cancel because of the equilibrium condition fi^ + fJ-p — fJ'N = 0. The substitution also 
causes — )• —m*^ such that we are left with a positive transition rate, 

PQ=^-^m*j^m*pfieT^I (4.20) 
where I is an analytic integral, 

1= r dxi r dx2 r dxs — (^1+^2+3:3)^ — ^ 

Loo J-ix,+.,) =^(l + e-i)(l + e-2)(l + e-3) ^ ' 

= I (vr2C(3) + 15C(5)) (4.22) 
« 20.56. (4.23) 

Putting everything together, and using the numerical estimates from the nuclear matter 
discussion, the transition rate becomes, 

w ^Gl 
n 87r5 



^ - (1 + ) m% ml /ie / (4.24) 



3-5 X 10^^ ( -^^) iT,f cm-3 (4.25) 
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where Tg = r/(10^ K) is the dimensionless, scaled temperature. A typical value for the 
reduced mass factor is 0.8. The temperature dependence is consistent with literature — 
remember that we are calculating the transition rate, not the luminosity. The luminosity 
contains an additional factor of energy in the integral that contributes an extra factor of 



4.1 Estimate of the current from direct Urea 

The first step is to ensure that the electron can actually escape the star before it decays via 
the weak interaction. This involves calculating the mean free path and seeing if it is larger 
than the radius of the star. In the literature for calculations of this type, such as neutrino 
luminosity for cooling, $7 is the volume of the neutron star — it is necessary to account for 
all the transitions that occur in the entire star. Here, we are interested in the decay rate 
of a single electron, so we will take O to be the volume in which a single electron exists, 
which is the inverse of the electron number density. 



Assuming that the electron travels at the speed of light though the protons — due to 
the non-dissipative nature of the current with respect to the electromagnetic interactions, 
the mean free path can be found using, 



The typical radius for the neutron star is i? ~ 10 km. We see that for T < 10^" K the 
electrons can easily escape the degenerate region before the P-odd asymmetry gets washed 
out due to the weak interactions. The counterintuitive density dependence — at higher 
densities the mean free path is larger — is a natural consequence of Pauli blocking. As the 
density of the star increases the number of protons increases and we would expect a shorter 
mean free path, but the number of neutrons increases as well. The suppression due to the 
higher neutron chemical potential is greater than the enhancement gained by increasing 
the proton density. 

Using (3.3) we are now in a position to estimate the magnitude of the current travelling 
along a single quantum of flux. 



There is no dependence on the magnetic field because we have normalized per unit of 
quantum flux. As discussed earlier this reaction is dominant when hyperons appear at 
n/iiQ > 3 or when the proton fraction is large Xp > 1/9. 

From (4.27) we see that is it much easier for an electron to escape a cool star rather 
than a hot, newly born star, T > 10^^ K. If the star is very hot the electron is not able to 
keep its asymmetry due to the weak rescattering, as one can see from (4.27). As it cools 



r. 




(4.26) 




(4.27) 




(4.28) 
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there will be a critical value where the electrons can escape, but are still created at a large 
rate, meaning the current is very large (4.28). The current is largest when the star is hot, 
but not so hot that the electrons are unable to escape the region of degeneracy, fie ^ T. 
This temperature is roughly determined by the condition that the electron mean free path 
with respect to weak interactions is approximately equal to radius for the neutron star is 
i? ~ 10 km. It is expected that this temperature drastically depends on the equation of 
state and other specific properties of the environment as (4.27) suggests. 

4.2 The effect of a large magnetic field on the transition rate 

In this short subsection we will argue that the effects of the magnetic field can be safely 
neglected in calculating the transition rate. Of course, the magnetic field still plays a crucial 
role in producing the required asymmetry for the current, see Section 3 and the Appendix 
for details. There has been much consideration of Landau levels on the rates of processes 
that occur in neutron stars. Landau levels can have drastic effects on the spin of electrons, 
though these effects are suppressed by the unusually high chemical potentials found in 
neutron stars. An extremely large magnetic field is needed for these effects to manifest 
themselves in neutron stars. It must be roughly comparable to the chemical potential 
eB ~ /ig in order for substantial changes for the transition rates to occur. Numerically, 
this corresponds to very large fields B ~ Bcii^/ml ~ 10^''G, which are much larger than 
the magnetic fields found in typical neutron stars. 

To introduce Landau levels to the earlier calculation we can simply replace the electron 
dispersion relation with. 



where n are all natural numbers, but only electrons with a spin antiparallel to the magnetic 
field are allowed in the lowest, n = 0, level and h = B/B^ is the magnetic field normalized 
to the critical field, B^ = mg/e = 4.4 x 10^'^ G. The the electron phase space becomes. 



where the term after the sum is the degeneracy per unit area and gn, equal to 1 for n = 
and 2 otherwise, counts the spin degeneracy. The maximum Landau level occurs for pz = 
when all the energy goes to putting the electron in the Landau level and none goes into 
the momentum. This yields. 



If the number of Landau levels is very large, nmax ~ 2bm^ ^ which is a common case for 
a typical neutron star, then the phase space returns to the the one we used in the previous 
section (4.24). Therefore, for a typical neutron star the transition rate basically remains 
the same as the electron phase space essentially unchanged, 



El=pl+ml{l + 2nb) 



(4.29) 




(4.30) 






(4.32) 
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where we used Ef. ~ /ig and ^Ug ^ nie- This result is in accordance with our previous rough 
argument that a very large magnetic field is required to produce any substantial changes. 

We want to contrast this generic case with a rare situation when only the lowest Landau 
level is accessible, rimax = 0. This occurs when the neutron star is abnormally cold, the 
magnetic field is unusually large, or the electron chemical potential is unnaturally small 
[41]. In this approximation we recover the usual electron dispersion relation, while the 
available phase volume becomes, 

n I 0^^mlt. (4.33) 

As before, we move the momentum dependent parts of the phase space outside of the 
integral by approximating them by their Fermi momenta. The rest of the integral takes 
place identically. We can compensate for the magnetic field in our calculation of the 
transition rate simply by taking, 

b 

fie^ — - (4.34) 

In fact this transformation can be done to account for the magnetic field for any of the 
other rates we derive (5.19), (6.9), (7.7). In particular, for the direct Urea case we are left 
with, 

^ = g(l + 3Cj)™^»p!^T=/. (4.35) 

which we can compare to equation (4.24). This expression is valid in three cases: T <^ 10^ 
K, He ^ 100 MeV, or 5 ^> 1. Only the last of these has a real chance of manifesting itself; 
magnetars have huge magnetic fields B ~ 10^^ G. Otherwise the Landau levels do not have 
a significant contribution and the earlier expression is valid. 



5. Nuclear Matter: The Modified Urea Process 

The direct Urea process is the simplest to consider, but it is unlikely to be the most 
common process in a star. It is much easier to conserve momentum if a nucleon is included 
to supplement momentum transfer. We will consider the mean free path of an electron 
scattering of a proton assisted by a neutron, equation (2.10). The inclusion of the nucleon- 
nucleon interaction into the matrix element is non-trivial and we will use the one pion 
exchange/landau liquid method found in [42], the interaction illustrated in Figure 4. 

There are many other diagrams similar to the one illustrated. They include all pos- 
sibilities of pion exchange between protons and nucleons as well as crossing diagrams to 
account for the final neutrons being indistinguishable. Summing them all causes the vector 
contributions to cancel and we are left with the scattering matrix approximated to be, 

Ei^r^-^Sci (£)*„„„., ,5.1, 
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where / ~ 1 is the p-wave vrA^ couphng constant, Ca_ = 1-26 is the Gamow-TeUer couphng, 
and aurca ~ (0.63 — 1.76)(no/n)^/'^ is a factor that accounts for the pion propagator and 
the short range Landau hquid contributions. Fohowing [42], but somewhat preempting the 
calculation, we approximate the propagator of the internal nucleon using the total lepton 
energy /ig — remember that the energy the neutrino carries away is negligible. 
To calculate the mean free path we will first calculate the transition rate. 



w 



n I d^Pi S {27T)'n S^^\pf - Pi) J2 \^\^ ' (5-2) 



(2vr)i8 , 

i=l spins 



where S is the Pauli blocking factor. 

We separate the transition rate into angular and radial parts leaving. 
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(5.3) 



(2vr) 

spins 



where, 

6 

P = l[jpUPiSdiEf-Ei) (5.4) 

i=l 
6 

Q = Y[J dni6'^^\pf-p^^. (5.5) 



i=l 



We start with the Q integral. The contribution from the momentum of the neutrino is 
negligible compared to the rest so we neglect it in the (^-function. This allows us to take 
an integral over all angles causing the term with angular dependence to vanish. We then 
follow [34] in doing the angular integrals to get, 

Q = 7^ • (5-6 



-27- 



We can now do the PQ integral. Changing variables from momentum to energy and 
approximating the energy by the particles fermi energy changes the measures to, 



Pidpi 


= m*^dEi 


pldp2 


= fie m*pdE2 


Pldps 


= l^ldE^ 


Pidpi 


= m*j^dE4 


vldpb 


= m*j^dE^ 


pldpe 


= EgdEQ , 



(5.7) 
(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 

where we took kp ~ /ie- Then we take the Eq integral over the delta function and make 
substitutions Xi = ziz/3{Ei — /ij) to make the remaining integral dimensionless as in the 
direct Urea case. Compared to direct Urea, the two extra particles contribute two more 
Fermi distributions to the statistical factor, 

1=1 ' 

Using the beta equilibrium condition, fip + fie = all chemical potentials cancel, leaving, 

PQ = M! ^rn*pf (m^)4 fi^ I (5.14) 

where I is the integral, 

5 

r 



TT^^) e(xi + X2 + X3 + X4 + Xs) (5.15) 

4 / i-oo \ poo I ^ \ 1 

=n / ^-0/, Ex, nrr^ (^-^6) 

I— I \J~oo / J -(XI+X2+X3+X4,) \j=l I k=l ^ 

= ^C(7) + f ^^C(5) + ^vr^C(3) (5.17) 
«192. (5.18) 

Gathering all the terms together we get the transition rate per unit volume for an electron 
scattering off a proton assisted by a neutron, 

n = —9 ^4 «U.ca T I (5.19) 

= 9.2 X 10^6 ( f'e \ (!nk\\7 ^^-1 ^^.3 (5 20) 

VlOO MeV/ \mN J 

where as discussed earlier ~ m*p ~ 0.8mAr. The rate of this process is much smaller 
than we found earlier in the direct case (4.24) and the temperature dependence is now 
rather than T^. Also, because of the density dependence of cturca? the transition rate is 
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independent of density. Both of these differences result from this being a higher order cal- 
culation, which involves two extra particles — the ingoing and outgoing neutron assisting in 
momentum conservation; the transition rate is suppressed by a factor of the QED coupling 
constant, and the two extra factors of T come from the two extra measures of integration. 

5.1 Estimate of the current from modified beta decay 

We are interested in the mean free path of an electron. The volume occupied by a single 
electron is Qe ~ 10~^^ cm^, and the electron moves at the speed of light, leaving the mean 
free path to be, 



At the beginning of the star's life, when it is very hot T ~ 10^^ K, the electrons are 
trapped, then as the star cools the electrons are allowed to escape. Here we get a similar 
counterintuitive density dependance as (4.27) but the effect is stronger because there are 
two neutron Pauli blocking terms to contend with. 

We now follow the same prescription as earlier to derive the current. 



This is by far the weakest current from any phase of matter. As before the the current is 
strongest early in the star's life when it is hot, and though it is suppressed it can get quite 
large at high temperatures due to the dependence. 

6. Kaon Condensate 

As the density of the star gets above three times nuclear density is it possible that a charged 
kaon condensate will appear [35] . It is now energetically favourable for electrons to scatter 
off the condensate and turn into neutrinos. We are interested in calculating the transition 
rate of an electron decaying in the presence of a kaon condensate, equation (2.12). The 
effect of condensates on the scattering matrix was first used to describe pion condensates 
[37] and later for kaon condensates [39, 43] and involves evaluating the process given by 
the Feynman diagram in Figure 5. 




(5.21) 




(5.22) 




Figure 5: Electron decay in a kaon condensate. 
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The process needs the aid of a neutron quasiparticle, similar to that of the modified 
beta decay, but not because it is required to conserve momentum. We start with the usual 
hadron and lepton currents and the kaon condensate is described by a chiral rotation in 
V-spin space, U = e^'^^^ . For kaons a small rotation can have a large effect so only 
must be considered. The matrix element, as found in [39], is given by, 

El'^P = ^(l + 3Ci)sin2 0,, (6.1) 

spins 

where ~ 0.1 is the kaon amplitude, Ca = 1-26 is the Gamow- Teller coupling, and 
9c ~ 13° is the Cabibbo angle. 

We are interested in the regime where the kaon momentum p3 is zero. The transition 
rate is given by. 



w 



(2vr) 



12 



jd^pid^P2d^P4d^P5S{27r)'^n5{pf - pi)\S\'^ . (6.2) 



where S is the statistical factor containing Fermi blocking terms, is the volume of the 
neutron star. The subscripts on pi, p2, Pa, and ps, label the momentum of the ingoing 
neutron, the electron, the outgoing neutron, and the neutrino, respectively. Following the 
work done in [34] we can separate the angular and radial integrals so each can be evaluated 
independently, 

q5 

^=(^El'5|'^Q> (6-3) 

spins 

where, 

P = j p\dpi pldp2 pldpApldp5 S6{Ef - Ei) , (6.4) 

Q = j dnida2d^idn^5^^\pf -Pi) . (6.5) 
We start by doing the angular integral Q is done the same as in earlier cases, 

2piP2P4, 

The PQ integral can now be started, 

PQ = TT^^ / Pldpi pldp2 pldpi P5dp5 S 5{Ei + E2 - E3 - E^ - E^) . (6.7) 

We first change the variables of integration from momentum to energy pidpi = EidEi and 
perform the neutrino integral over the delta function. In an effort to make the final integral 
tractable, we follow [34] in approximating factors next to the measures as constant. The 
neutrons are non-relativistic so their energy is just their effect mass, m^, and the electron is 
ultra relativistic, because of its large Fermi momentum, and the energy is just its chemical 
potential, /ie- As before, these factors can then be moved outside the integral. It is also 
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convenient to change variables to facilitate the final integral over the Pauli blocking factors. 
Changing variables takes {Ei + E2 — E3 — E^Y to {xi + X2 + X4)^r^, where we used the 
equilibrium condition /i^ = /ie- Upon changing variables the measure of integration gives 
us a factor T'^ and we are left with, 

PQ=^^{myfi,,T'l (6.8) 

where I is the same analytic integral from the direct Urea process (4.21). 

Putting everything together we get the transition rate per unit volume of an electron 
decaying into a kaon, 

^ = ^(l + 3Cl)sin2 0,("^^)VeT5/ (6.9) 

= 4.4 X 10^9 ( \ f!I^y f -V^'(r9)^ s^^ cm-3 (6.10) 

\imMe\ J \mN J \noJ ^ ^ 

where Tg = T/(10^ K) is the scaled temperature, and ~ O.SmTv- The temperature 
dependance is the same as the direct Urea process (4.24), but we would expect it to be of 
higher order than direct Urea because it involves an extra particle, and thus smaller. The 
temperature dependence is the same because of the way we treat the kaons as a condensate, 
rather than an extra particle. The condensate picture is a rotation of the direct process, 
rather than a whole new particle interaction as in the modified Urea process. 

6.1 Estimate of the current in a kaon condensate 

Because of the high chemical potential the electrons must be relativistic leaving the mean 
free path to be, 

/ \ 4/3 

4 = 3.0 X 10^2 (j.g)-5 / ^ \ (g_^-L) 

Once again the mean free path larger than the radius of the neutron star and the electrons 
can escape. We also notice the counterintuitive, but now familiar, density dependence 
discussed in the direct Urea case. This effective mean free path with the helicity intrinsic 
in the weak interaction creates a current given by equation (3.3), 

(j) ~ 3.6 X 10^" (Tgf (^j MeV. (6.12) 

Though the generation of the current in a kaon condensate happens though a wildly differ- 
ent process than direct Urea, and the numbers we use are quite different, we see that after 
the star has cooled the numbers conspire to give currents of similar magnitude (4.28). 

7. Quark Matter 

The last case we consider is what would occur if the hadrons separated into their constituent 
quarks. This is the case in quarks stars, where degenerate quarks exist deconfined. In 
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calculating the mean free path of electrons in quark matter there are two possible reactions 
we must consider, if we restrict ourselves to the lightest quarks, given in equation (2.13). 
Each is given by the Feynman diagram in Figure 6, where the d quark can be substituted 
out for the s quark. 




Figure 6: The direct Urea process for quarks. 



Unlike nuclear beta decay, the lowest order quark beta processes do not require help 
from an external particle to to conserve momentum; they proceed unsuppressed [36]. Be- 
cause they are deconfined, the Fermi momentum of the quarks are much closer to each 
other than the Fermi surfaces of the hadrons in neutron stars. We will first consider the 
transition into the down quark. The transition rate is given by, 
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d^pi d^p2 d^p3 d^pi S {2Tifn {pf-pi)\S\ 



(7.1) 



where S is the statistical factor containing Fermi blocking terms, is the volume of the 
phase space. The subscripts on pi, p2, ps, and p4, label the momentum of the ingoing 
electron, the up quark, the down quark, and the neutrino, respectively. The matrix element 
of this process is found in [36], 



^G^cos e^— 



1 



P3 P4 



(7.2) 



which looks similar to the usual four Fermi matrix element but contains correction due 
to the ability to relate the quark's momentum to the quark gluon coupling constant by 

2 

= 1^. We have included the factor to account for the degrees of freedom of the quark 
in here. The transition rate can be split into two integrals. 
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(7.3) 



(7.4) 



(7.5) 



These integrals are nearly identical to the kaon case. The Q integral is the same as in the 
direct Urea case (4.9). We make an approximation to the remaining PQ integral when we 
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change variables from momentum to energy. As we have done in previous integrals, the 
the momentum of the electron and quarks is replaced by their value at the fermi surface. 
Performing the neutrino integral, di?4, over the delta function leaves us with the same 
integral as equation (4.15). We perform a similar change of variables as before to make the 
integral dimensionless. Following the same procedure as the direct Urea case, the quark 
beta equilibrium condition /ig + /^m = fJ^d causes the chemical potentials to cancel, and the 
PQ integral becomes, 

PQ=^-^fiekukdT'>I, (7.6) 

where I is the same integral given in equation (4.21), which is the same as the direct 
Urea and kaon cases. An identical calculation can be done for the electron scattering into 
strange quarks. Putting everything together we get the transition rate for electrons in 
quark matter, 

Wd = 0,^Gp cos^ Ocas kekukdT^ I (7.7) 




To estimate the the magnitude of the transition we assume that kd ~ ku ^ kg as in 
the massless noninteracting case and use the relationship ^ fi^ = {3Ye)^/^kq, where 
Ye = ne/rih is the ratio of electrons to baryons in the star. The Fermi momentum of the 
quarks is estimated using (2.23). We follow [36] in estimating ~ 0.4 and Ye ~ 0.01 
The total transition rate is the sum of Wg and Wd, 

]^ = 1.8 X 10^1 ( ^] {Tgf cm-3 (79) 

We see the temperature behaviour that has become a signature for the first order, four 
Fermi interaction. 



7.1 Estimate of the current in quark matter 

As before, the first step in finding the current is estimating the mean free path of the 
electon in quark matter by assuming that the electron propagates at the speed of light. 
The same physics creates the current in quark matter and it is still non-dissipating. We 
now use the new definition of Ue to determine the volume occupied by a single electron, 
17 = {Ye ^^^,)~^, and we are left with, 

4 = - = 2.9 X 10^° (rg)-^ km (7.10) 
w 

The radius of the star, i? ~ 10 km, is much smaller than this and the current will propagate. 
In quark matter the mean free path is not dependent on density; the Pauli suppression and 
enhancement cancel each other — see section 4.1 for a brief discussion. When the electrons 
reach the crust they are removed from the system creating a new effective mean free path 
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for the electron. This effective mean free path with the heUcity intrinsic in the weak 
interaction creates a current, 



The typical density for quark matter is nb ~ 10 no, but could easily be higher. Once again 
the numbers have conspired and the magnitude of the current is close to the value for 
both direct Urea (4.28) and kaon (6.12) processes. They are all first order processes, but 
we see a critical difference in the density dependence. Unlike the other currents the quark 
current actual gets larger with increasing density. This happens because quark stars remain 
charge neutral in a fundamentally different way than neutron stars (see Section 2.4) and 
the electron chemical potential is determined differently. 

8. Applications in Neutron Stars 

In this section we discuss five different applications of the topological vector current intro- 
duced in this paper. 



Section 8.1: Neutron star kicks 
Section 8.2: Pulsar jets 
Section 8.3: Toroidal magnetic field 
Section 8.4: Magnetic helicity 

Section 8.5: Type-II superconductivity and precession 



These phenomena are observed in many neutron stars and appear to be unrelated to 
each other, but we will argue that they originate from the same source — non-dissipating 
topological currents. All these phenomena have one thing in common: they are P-odd 
effects. 

We have calculated the magnitude of the topological current in four types of nuclear 
matter and summarize them in Table 1. Direct beta decay, which is also the preferred 
reaction when hyperons are present, kaon condensates, and quark matter, all are radically 
different, but they all make a current about the same order. There is a stark difference 
in the modified Urea process, which occurs in ordinary nuclear matter where the proton 
fraction is below 1/9. The current is four orders of magnitude smaller than the rest. This 
is because it is a higher order process involving more particles. It results in a more severe 
dependance on temperature which suppresses the process when the star cools. 

The estimates of the current are relatively small in magnitude, but much larger currents 
exist on microscopic scales. The estimates represent a component of the current that 
produces a coherent effect in the entire region where matter is degenerate and the magnetic 
field all points in one direction. As a technical remark, the estimates of the current {j) 
are presented in MeV units. The current is defined as a number of particles crossing the 
surface equivalent to one quantum of magnetic flux per unit time. We can obtain the 
electromagnetic current in conventional units by multiplying by e = y/Ana and converting 
the result into Amperes: e{j) ~ 10^ A for {j) = 1 MeV. 




(7.11) 
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Table 1: A summary of the current per magnetic flux quantum calculated from each interaction. 
The value is for a single quantum of flux, such as that found in a type-II vortex. The appropriate 
density for the kaon and the direct Urea currents is n = 3 no, for quarks is nb = 10 no, and for 
modified Urea is n = ng. 



Current {j) (MeV) 



Direct Beta Decay 3.0 x 10"^ (^) (rg)^ 



4/3 



Modified Beta Decay 7.7 x lO^^^ (^^j (rg)^ 

/ \ -2/3 

Kaon Condensate 3.6 x 10"" ( (rg)^ 
Quark 2.8 X 10-9 (^)^^^ (rg)5 



8.1 Neutron star kicks 

We will briefly discuss the application of topological currents to neutrons star kicks. A 
much more detailed analysis can be found in [23]. It is accepted that pulsars have much 
higher proper velocities than their progenitors, some moving as quickly as 1000 km/s 
[24, 44, 45, 46, 47, 48, 49, 50]. Such large velocities are unambiguously confirmed with the 
model independent measurement of pulsar B1508+55 moving 1083l^[!^ km/s [22]. There 
also appears to be a strong correlation between the direction of the kick and the rotation 
axis of the star [26]. If the electrons carried by the current can transfer their momentum 
into space the current could slowly and steadily pushing the star over time, much like a 
rocket. This is contrary to other mechanisms where the kicks happens shortly after the 
star's birth. The sustained nature of the kicks allows a rotation-kick correlation {P • Q) to 
form regardless of the angle between the star's rotation and magnetic field. Neutron stars 
spin on the order of milliseconds and the kick occurs over hundreds of years. The force 
vectors from the kick form a cone aligned with the axis of rotation that average to produce 
a kick along the axis of rotation. Long duration kicks are supported by the analysis in [25]. 

Currently no mechanism exists that can reliably kick the star hard enough. Explosions 
during collapse can only reliably kick a star 50 km/s [24], asymmetric explosions can only 
reach 200 km/s [50], and asymmetric neutrino emission is plagued by the problem that 
at temperatures high enough to produce the kick the neutrino is trapped inside the star 
[51, 52]. More exotic mechanisms require exotic particles and more fine tuning. We pro- 
vide a rough calculation demonstrating that topological kicks can provide the momentum 
required and do so because of the immense electron chemical potential, not because of 
the temperature of the star. A detailed calculation of the kick is presented in [23] . Also, 
because of the correlation between the magnetic field and the spin axis, topological kicks 
naturally align with the spin axis of the star. 

Though it is unclear what will happen to the electrons when they reach the surface of 
the star, we can estimate the size of the kick if we assume the entire momentum carried by 
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the topological current will be transferred into space by some means. This assumption is 
likely correct for bare quark stars where the crust (the region where ~ T) is only about 
1000 fm wide [53]^ and is most likely wrong for typical neutron stars where the crust is 
about 1 km thick. The energy of the electrons would be absorbed by the crust and would 
not contribute to the kick. As a consequence we conjecture that stars with very large kicks, 
V ^ 200 km/s, are quark stars and that slow moving stars, v < 200 km/s, are kicked by 
some other means [50] and are typical neutron stars. Confirmation of this classification 
would provide a simple and well formulated principle that distinguishes quark stars from 
typical neutron stars and results in a bimodal distribution for kick velocities as supported 
in [49]. As is known, other criteria such as mass, size, cooling rate, etc. cannot easily 
discriminate between quark stars and neutron stars, see e.g. [33, 56] for review. 

The first step in estimating the magnitude of a topological kick is determining the total 
momentum transferred to the star. There are A'v ~ 7 • 10^^ B/Bc quantum units of flux in 
a star that can be distributed in either superconducting domains or vortices, see eq. (2.26). 
The current is independent of the internal structure of the star. If all the electrons are shot 
out of the star or transfer their entire momentum, then the momentum is given by the total 
current (3.1), the number of electrons that leave the star per unit time. The momentum 
a single electron transfers out of the star is equal to its Fermi momentum, ke = 100 MeV. 
The fuel for the kick is the chemical potential, not the temperature. Therefore the kick 
may continue even when star is already cool; the kick in our mechanism is not an instant 
event, but is rather a long, slow, steady process that pushes the star. Putting this all 
together, the current transfers A^v (i) units of momentum per unit of time. 

The appropriate way to estimate the magnitude of the kick is to integrate the rate of 
momentum transfer per unit time by taking into account the time evolution of the star as 
the current is very sensitive to the temperature/time as Table 1 demonstrates. The rate 
of momentum transfer also changes as the star cools and the environment (density, phase) 
changes. We neglect all these complications and take a value for the current corresponding 
to T 10^ K, which is far from the maximum current can potentially be induced. 

The momentum transfer required for the star to reach a velocity oi v = 1000 km/s 
is quite large. Per baryon, the momentum required is rUn v ~ 3 MeV. The total baryon 
number of a neutron star of one solar mass is B^ ~ 10^^. The momentum for the entire 
star is then P ~ 3 • 10^"^ MeV. If we choose a current (j) = 10 (Tg)^ MeV, corresponding 
to an average current in Table 1, the time required to attain this momentum is, 

N^keU) VlOOOkm/s; U A (i) ) \ T J ' 

w 10,000 years, (8.1) 

We restore the canonical dimensionality of MeV~^ in seconds by multiplying h = 6.58- 10~^^ 
MeV-s. This is a conservative estimate of the kick strength. The current is much larger at 

^It has been recently argued that a crust in quark stars could be much larger in size than previously 
thought due to development of a new heterogeneous mixed phase [54, 55]. As we mentioned above, it is not 
our goal to discuss the interaction of current with the crust, however, an intense current from the core of 
quark star may destroy the crust in this new mixed phase in few locations similar to volcanoes on earth. 
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the star's birth (when it is hot) and kicks such as those seen in, for example, the Vela pulsar 
can be easily explained. If electrons actually leave the star, rather than transferring their 
momentum through radiation, the electron chemical potential will slowly decrease and the 
current may stop running. Charge neutrality will cause matter to accrete isotropically and 
possibly maintain some of the chemical potential. 

As previously mentioned it is unlikely that the electrons could escape a typical neutron 
star — the crust has a thickness on the order of 1 km and electrons leaving the degenerate 
core would quickly be absorbed. On the other hand, quarks stars have a thin interface 
between the degenerate matter and open space [53] that electrons might be able to easily 
penetrate, or transfer their momentum by emitting photons into space. This difference in 
crusts could provide a valuable distinction: stars with large kicks are quark stars, stars 
with small kicks are neutron stars. 

Our kick mechanism is similar to neutrino driven kicks in that both use P-odd effects 
and particles to carry momentum out of the star. The electron and neutrino have simi- 
lar mean free paths, but the electrons leave much more slowly than the neutrinos. The 
fundamental difference in the two carriers. At low temperatures, when the neutrinos can 
escape the star, they do not carry enough energy to explain the kick. In contrast, the elec- 
trons that make the non-dissipating current (or more precisely the quasi-particles which 
freely travel along the magnetic field) carry very large momentum ~ /Xg. As a result the 
neutrino carries too little momentum when the mean free path becomes sufficiently long, 
while the momentum carried by the topological current remains very high even at very low 
temperatures T ~ 10^ K. This is the advantage of our kick mechanism. 

8.2 Pulsar jets 

A different but likely related phenomena is the recent observation of pulsar jets [57] that 
are apparently related to neutron star kicks [58, 59]. It has been argued that spin axes and 
proper motion directions of the Crab and Vela pulsars are aligned. Such a correlation would 
follow naturally if we suppose that the kick is caused by a non-dissipating current as we 
mentioned above. The current, and thus the proper motion, is aligned with the magnetic 
field, which itself is correlated with the axis of rotation. It would be very tempting to 
identify the observed inner jets [57] with the electrons/photons emitted as a result of the 
induced current. In this sense the mechanism for the kick is similar to the electromagnetic 
rocket effect suggested previously [59]. 

It is possible that evidence of the topological current may be directly detected in 
these jets. An observational consequence of the current is that a component of the X-ray 
emission in the trail of the neutron star will be left circularly polarized. Because of the 
parity violation in the star, only left-handed electrons will contribute to the kick. When 
these left-handed electrons interact they will create mostly left-handed photons, which 
coherently will be seen as left circularly polarized X-rays. As there are many other sources 
of X-rays the contribution is likely very small, but it further motivates the need for higher 
precision X-ray polarimetry [60]. 
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8.3 Toroidal magnetic fields 

There is a strong theoretical evidence for the existence of toroidal fields in neutron stars 
based on the stability of the poloidal magnetic field. References [29, 28, 30, 61, 62] argue 
that toroidal and poloidal fields of similar magnitudes must be necessary to stave off hy- 
drodynamic instabilities — the toroidal field suppresses poloidal instabilities and vice versa. 
Much has been done to find the observational consequences of a toroidal field, for example 



Estimating the toroidal magnetic field is a very complicated problem that requires 
a self consistent solution of the equation of the magnetic hydrodynamics. Our induced, 
topological currents represent only a small part of the system. We are not attempting to 
solve this problem. Instead, we shall argue that the currents we estimated are more than 
sufficient to induce the toroidal magnetic field correlated on large scale of order 10 km. 

A natural consequence of having a current running parallel to the poloidal magnetic 
field is that a toroidal component i^tor will be induced. The size of the field can be 
calculated naively using Ampere's law, but there is a subtlety because the magnetic field 
is being induced inside a superconductor. The magnetic field observed in neutron stars 
B ~ 10^^ G is actually induced by a much larger field H. The suppression comes from 
the perfect diamagnetism of the proton superconductor (the Meissner effect). This perfect 
diamagnetism is ruined at a critical field He ~ {^o/ AttX^) ~ 10^^ G where flux penetrates 
the star through small regions where superconductivity has been destroyed (vortices or 
domains). The supercurrents responsible for the perfect diamagnetism do not flow as 
easily and a small field is induced. 

Regardless whether the flux penetrates the superconductor as single vortices or in flux 
domains, we can assume that the superconductor is type-II.® The relationship between the 
applied magnetic field H and induced magnetic field S in a type-II superconductor is very 
nonlinear. The details have been worked out in [64] and [65], where the latter is a direct 
application to neutron stars. The important points are that below the first critical field 
H < Hci there is no magnetic field B induced. Just above the critical field a magnetic field 
appears that is approximately B ~ lO^^/Zci- As the applied magnetic field is increased 
above Hd the induced field starts to approach the applied field. 

We want to determine if the topological current produced by the poloidal field can 
induce a sufficient toroidal field by finding the length scale where Htor ~ H. Following [65] 
we assume that H ~ and we get the relationship B ~ \Q~^H for our magnetic field. 
We apply Ampere's law for a region of size L, 



V "^0 / 

where the expression in brackets describes the number of unit fiuxes bundled in the area 
TT-L^ such that we get the total current enclosed in our loop. We take $0 = '^/e and 

*The mechanism for type-I like superconductivity discussed in [10] relies on the electromagnetic inter- 
action between currents carrying vortices, not in altering the value of the Landau- Ginzburg parameter 
K = \/£,. We then still use results from type-II superconductors (indicated by k > l/v^) but the vortices 
are now bunched together in large domains with higher winding numbers. 



[63]. 
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substitute use our relationship between H and B. 

The naive estimate leads to the following expression for i^tor in terms of magnitude of 
poloidal magnetic field H, 

"'-~Mi)L~4^Ji^M^)- (8.3) 



This shows that a typical current from Table 1 can induce a toroidal field the same magni- 
tude as the poloidal field on scales the order L ~ 1 km, within the typical size of a neutron 
star. It is quite obvious that our estimate becomes unreliable when -fftor ^ H and we 
can no longer ignore the current induced by the toroidal field. For i^tor ^ H the problem 
requires a self consistent analysis which is beyond scope of the present paper. The point is 
that the toroidal field obviously develops as a result of topological currents and eq. (8.3) 
shows that its magnitude can easily become the same order as the poloidal magnetic field. 
If superconductivity is completely destroyed, as we will discuss Section 8.5, B = H and 
the toroidal field can induced on a much smaller scale, L ~ 1 m. 

8.4 Magnetic helicity 

Magnetic helicity, see e.g. [66], is defined as, 

■H = jd^xA-B. (8.4) 

The magnetic helicity is a topological object that can be expressed in terms of the linking 
number ,^2) of t'^o curves 71 and 72. The precise relation between T-l and interlinked 
flux $1 and <I>2 is given by, 

n = 2^>i$2 = N1N2 (8.5) 

where $1 = ^I'o-^i and ^2 = ^0-^2 are expressed in terms of unit flux $0 and the linking 
number is simply reduced to ?T'(7i,72) = N1N2. Therefore Ti takes integer values up to a 
normalization 2<I>g. This linking number is preserved, ^ = 0, in a magneto- fluid with 
zero resistivity, which is a very good approximation for neutron stars. This topological 
invariance provides the stability necessary for the poloidal field. 

We want to emphasize that the magnetic helicity is the dot product of a vector and 
a pseudovector, making it a pseudoscalar. Under the parity transformation x — t- —x the 
magnetic helicity is P-odd: Ti — )• —Ti. This implies that the magnetic helicity can be 
only induced if there are parity violating processes producing a large coherent effect on 
macroscopic scales. Many attempts to generate helicity rely on instabilities in the magnetic 
field caused by the star's rotation. Such correlations B ■ $1 are P-even, and though they 
may generate toroidal fields they cannot be responsible for helicity. 

Our observation here is that the non-dissipating topological current introduced in 
the present work has precisely this property: the topological current produces the P-odd 
correlation (P-B) and is capable of inducing magnetic helicity Ti ~ j. In fact, our estimate 
for the induced toroidal field /?tor unambiguously implies that the magnetic helicity will 
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be also induced, see eq. (8.3). The magnetic flux from the toroidal and poloidal fields is 
always interlinked and contributes to the magnetic helicity, 



where ^ and ^I'torr are the original poloidal and induced toroidal magnetic fluxes (8.3) 
correspondingly. 

Strong observational evidence, see [63] and references therein, supporting the presence 
of the toroidal component unambiguously suggests that the magnetic helicity Ti must be 
non-zero in neutron stars. The P-odd quality of the magnetic helicity may be strong, 
indirect evidence supporting our claim that P-odd topological currents have been induced 
at some moment in the star's life. Otherwise, it is very difficult to understand how such a 
large, coherent P-odd effect could be produced. 

8.5 The conflict between vortices and precession 

It has been observed that neutron stars precess [67, 68] and that the degree of precession 
conflicts with the commonly held belief that the protons form a type-II superconductor in 
the core [16, 17]. It has been shown that if the superconductor is type-I there is no conflict 



When a magnetic field is applied to a type-II superconductor the flux finds it energet- 
ically favourable to penetrate it by forming many vortices each carrying a unit of quantum 
flux. In a neutron star the large number of these get tangled with the superfluid neutron 
vortices that have formed to carry angular momentum. If the star precesses with a large 
enough angle the superfluid vortices must break through the superconducting vortices for 
rotation to continue. Incredibly large amounts of energy are dissipated in this process 
which would cause the star to stop rotating. We conclude then that the superconductor 
must be type-I where the flux bunches in large groups organizing macroscopically large 
domains and there is room for the neutron vortices to move around. 

The problem is that the Landau- Ginzburg parameter for a typical neutron star indi- 
cates the superconductor is type-II. A solution to this is suggested in [10]. If a sufficiently 
large current runs along type-II vortices an attractive force arises that causes the vortices 
to bundle together like they would in a type-I superconductor even though the Landau- 
Ginzburg parameters indicate type-II behaviour. 

Even if the current is not strong enough to make the vortices attract each other it has 
also been argued [10] that the mere presence of an induced, longitudinal current, arbitrarily 
small, would destroy superconductivity, thus resolving the problem. In many condensed 
matter systems such kind of instability has been experimentally tested, see [10] for relevant 
references on test of this instability in condensed matter systems. This instability can be 
delayed for small currents or even stabilized due to impurities. But the lesson from these 
condensed matter systems is that when a current aligns with the magnetic field creating 
the vortex the properties of the vortex lattice are completely changed or destroyed. 

We expect similar behaviour in regions of the neutron star where both the Landau- 
Ginzburg parameter suggests type-II behaviour and longitudinal currents are induced. 




[20]. 
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While many features of the system are still to be explored, the point is that that even 
small topological currents along the magnetic field will likely destroy the vortex lattice by 
replacing it with a new, unknown structure, similar to the condensed matter experiments 
mentioned above. The exact state is not essential at the moment, only that the Abrikosov 
lattice is destroyed by longitudinal currents and the conflict formulated in [16, 17] is re- 
solved. 

9. Conclusion 

The goal of this paper is to argue that a persistent, topological current is induced in neutron 
stars. All the requirements are present: a large degeneracy /.tg S> T, an approximately 
chiral, Dirac-like spectrum at /Ug ^ rrie, and a large magnetic field B. It is an unusual 
but beautiful pure, quantum phenomena that has no analogue in classical physics. It is 
fundamentally new, which raises a question: have similar phenomena been studied? There 
are systems with the potential to manifest anomalous currents. More so, effects similar to 
those discussed in this paper have been experimentally tested in some condensed matter 
systems and are going to be tested at RHIC and there is a possibility that these effects can 
be experimentally tested in terrestrial laboratories. 

In condensed matter laboratories low temperatures and strong magnetic fields present 
no technical difficulties. The key is finding a system of quasiparticles with a Dirac-like 
spectrum. There are such systems: superfluid He^, high Tc superconductors with d-wave 
pairing, and graphene. Remarkably, in superfluid He^ the current analogous to our anoma- 
lous current has been observed, see reviews [5] and [7]. 

The relativistic heavy ion collider (RHIC) at Brookhaven also has potential. An ana- 
logue to the anomalous current has been used to predict a charge separation effect [11] 
and preliminary experimental results are supporting it. In this analogue each requirement 
for the current to exist in the neutron star has its complement: the role of the coherent 
magnetic field is played by the angular momentum L, which occurs at non-central nuclei 
collisions, and the the role of parity violating effects is played by the induced 6 vacua. The 
observation of the charge separation indirectly supports our prediction of induced anoma- 
lous currents. We should remark here that very strong interactions in each nuclei-nuclei 
collision event cannot wash out the produced asymmetry. This parallels our argument in 
Section 3 that strong electromagnetic interactions do not wash out the P-odd produced 
asymmetry and the relevant scale of the problem is the mean free path of the electron due 
to the weak interactions that are capable of washing out P-odd effects. 

Confirmation of our claim that anomalous currents exist in neutron stars would have 
an enormous effect on the physics of neutron stars. In particular, it may explain neutron 
star kicks and pulsar jets (see Sections 8.1 and 8.2) and give a way to discriminate between 
neutron stars and quark stars, it may shed light on the nature of the toroidal magnetic field 
required for stability of the poloidal field (see Section 8.3), and it may resolve the confiict 
between the observed precession of a neutron star and type-II superconductivity commonly 
believed to exist in the core (see Section 8.5). Topological currents also provide a source of 
finite magnetic helicity, a P-odd topological invariant that does not decay in a neutron star 
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environment. This may shed some hght on the origin of the strong, self-supporting system 
of toroidal and poloidal magnetic fields in neutron stars (see Section 8.4). This necessity 
of finite magnetic helicity is strong, indirect evidence that non-dissipating currents move 



In Section 8 we mentioned many apparently unrelated observational effects: neutron 
star kicks, toroidal fields, and magnetic helicity. These all have a P-odd symmetry and it is 
likely that they all originate from the same P-odd physics. The topological vector current 
introduced in this paper occurs because of parity violating effects. This current may be 
responsible for all of these P-odd phenomena. 
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A. Calculation of the Helicity A(/i, B, T) 

In order to calculate the average helicity of the electrons, we must consider the creation of 
an electron through nuclear beta decay in a large magnetic field and chemical potential. 
There have been many calculations of transition rates in magnetic fields, notably [69, 70], 
and more recently for both large magnetic fields and chemical potentials [41] and references 
within. We are looking for a specific property of the weak interaction, the helicity of 
electrons produced in a large magnetic field and chemical potential. The helicity of a 
particle is given by A = cr ■ p/|p|, where cr are the Pauli matrices and p is the particle's 
spatial momentum. The two eigenstates of the helicity operator correspond to the values 



where ^ is the rest spin vector of the particle. The expectation value of the helicity can be 
calculated by looking at the ratio of decay rates T with different helicity. 



The small scale of the interactions and the immense magnitude of the magnetic field 
make it necessary to consider the Landau levels the electrons decay into. Detailed analysis 
of the effect of the Landau levels on interactions while considering the complete electron 
wave function has concluded that the only significant change is in the energy of the electron. 



along B. 



P^ = ±|P 
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(A.3) 
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where b = B / is the ratio of the magnetic field and the critical magnetic field, B^ 



9 '^ 
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he 

The magnetic field also affects the phase space of the electron. It can only carry linear 
momentum in the direction of the magnetic field. The magnetic field also causes degeneracy 
of levels. With this taken in to account the phase space becomes, 



I 



(A.4) 

2-K 2-K ^ ' 

We will confine ourselves to calculating the helicity for an electron created through 
beta decay, one half of the direct Urea processes. Figure 3. We sum over everything but 
the electron spin, 

Y,\M\'' ^IQGlE^Ep [{l + 2,C\)E,{Ee-^e-U (A.5) 

-(1 - ■ (Pg - meSe)] , 

where Ca_ — 1-26 is the axial current constant, Gp is the Fermi coupling constant, 

m[m + E) 

and ^ is the unit polarization vector. 

In order to find the decay rate we use the modified electron phase space. To find the 
total decay rate it is necessary to sum over the probability of an electron appearing in the 
each of the Landau levels. The sum truncates where the energy of a Landau level is high 
enough to be disallowed by conservation of energy. Specifically we let nmax be the largest 
n that satisfies E'^ > m1{l + 2nb). The total decay rate is the sum of the decay rate of 
each level, 

r=^r„, (A.7) 

n=0 

where 

n^^^ = {El-ml)/ml2b. (A.8) 

The highest Landau level is reached when all of the energy goes into putting the electron 
in the highest landau level and none into the momentum p^- We see that increasing the 
magnetic field decreases the number of levels to which the electron has access. When 
b > £'g/2mg — 1/2 the electron only has access to the lowest Landau level meaning that 
all the electrons are spin down. 

We must also account for the non-zero chemical potentials of the protons, electrons, 
and neutrons. Because the electron and proton have large chemical potentials their phase 
spaces are constricted. The presence of a Fermi surface means that only higher energy 
electrons and protons can be created. This is modelled by multiplying the usual phase 
space by Fermi blocking terms 1 — f{E), where f{E) is the Fermi distribution. With these 
considerations the decay rate into a single landau level n is, 

dTn = — ^"^^!!^^'f^|M|2(27r)454( _ ^ _ ^ ) 

2E,^2Ee2TT 27r 2Eui2TTf 2Ep{2tt)-^^ ' ^ ^ ^ ^ ■'^^ ' 

(A.9) 
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There are a few quantities that naturally align themselves with the magnetic field. Firstly, 
all spins are either aligned or anti-aligned so if we choose B = (0,0, B^), then 



4 = (0,0,^e)and ^^ = (0,0,^;v). 



(A.IO) 



Also, because of the Landau levels, the electrons have linear momentum only in the direc- 
tion of the magnetic field, 

Pe = (0,0,pe)- (A.ll) 



The matrix element is reduced to 



J2 \Mf ^ SGlE^Ep [(1 + iCl)E,{Ee - Pe^e 
-{I -CA)\Pi.\{Pe-meSe) cos 9] , 



(A.12) 



where is the angle between the z-axis and the direction of the neutrino momentum. 
The integrals up to the final electron integral are straight forward. We are left with. 



dTn = Agn dpeiE^ - Ep - Ee 



1 - 



Pe^e 

E, 



(l-/e), 



where. 



(A.13) 



(A. 14) 



The term Pe^e = MPe\ gives us the helicity eigenstates. We can get the combinations 
required for finding the average helicity, 



rn(+)-r„(-) = 



rpo 

4Agn / dpe{E^ 
Jo 

po 



Ep — Ef 



^2\Pe 

' E, 



(l-/e 



dpeiE^-Ep-EeYil-fe) 



(A.15) 
(A.16) 







In order to do the final integral over the Fermi distribution we appeal to the Sommerfeld 
expansion, 



h{E)il-f{E))dE = h{E)dE - — — h{E) 



(A.17) 



E=^l 



The particles in a neutron star are in chemical equilibrium. The maximum energy available 
for the electron is equal to its chemical potential, Eq = ji, leaving the integral part of the 
Sommerfeld expansion to vanish. Physically this makes sense because the equilibrium 
processes only occur thermally — the transition rate at zero temperature vanishes. Doing 
our integral has been reduced to taking a derivative. 



r„(+) - r„(-) = - E, - 

r„(+) + r„(-) = - '-^'^'^^(E^ -E,-E. ^^ 



Eedpe 



(A.18) 

(A.19) 
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where remember that E^ipe-, n). We can sum over each of these to get the total decay rate, 
then take the ratio to get the average hehcity. The important values are E^^ — Ep = /ie ~ fee 
and the sum goes up to rimax = k1/2m1h. The details after this are largely uninteresting 
and is computed by doing the sum numerically. As a check we find that the helicity at 
zero magnetic field is (A) = —1 as we expect. Over the range of fields we are interested 
B = 10^^ — 10^^ G the helicity is surprisingly constant. We arrive at, 

(A) = -0.84. (A.20) 

The helicity is close to —1, but not so close that it doesn't warrant a comment. With a 
large magnetic field the electron is forced into Landau levels. The single spin down state in 
the lowest level sometimes forces the electron into a right-handed configuration to conserve 
momentum. This occurs when the proton is created with the same spin as the initial 
neutron. In this case the electron and the neutrino must have then have opposite spins. 
If the electron is forced by the lowest Landau level to be spin down then the antineutrino 
must be spin up. Being right-handed, the antineutrino moves up in the direction of its 
spin. In order to conserve momentum the electron must move down, forming a right- 
handed configuration. The electron is bullied into being right-handed by the neutrino. 
The effect is not absolute though, and the intrinsic left-handedness of the weak interaction 
wins out. 

References 

[1] E. Wittcn, Superconducting Strings, Nucl. Phys. B249 (1985) 557-592. 

[2] A. J. Niemi and G. W. SemenofF, Fermion Number Fractionization in Quantum Field Theory, 
Phys. Kept. 135 (1986) 99. 

[3] D. Kharzeev, Parity violation in hot QCD: Why it can happen, and how to look for it, Phys. 
Lett. B633 (2006) 260-264, [hep-pli/0406125]. 

[4] D. T. Son and A. R. Zhitnitsky, Quantum anomalies in dense matter, Phys. Rev. D70 (2004) 
074018, [hep-ph/0405216]. 

[5] M. M. Salomaa and G. E. Volovik, Quantized vortices in superfluid heZ, Rev. Mod. Phys. 59 
(Jul, 1987) 533-613. 

[6] A. Y. Alekseev, V. V. Gheianov, and J. Frohlich, Universality of transport properties in 
equilibrium, the goldstone theorem, and chiral anomaly, Phys. Rev. Lett. 81 (Oct. 1998) 
3503-3506, [arXiv: cond-mat/9803346]. 

[7] G. E. Volovik, Superfluid analogies of cosmological phenomena. Physics Reports 351 (2001), 
no. 4 195 - 348. 

[8] M. A. Metlitski and A. R. Zhitnitsky, Anomalous axion interactions and topological currents 
in dense matter, Phys. Rev. D72 (2005) 045011, [hep-ph/0505072]. 

[9] M. A. Metlitski, Currents on superconducting strings at finite chemical potential and 
temperature, Phys. Lett. B612 (2005) 137-146, [hep-ph/0501144]. 

[10] J. Charbonneau and A. Zhitnitsky, A Novel Mechanism for Type-I Superconductivity in 
Neutron Stars, Phys. Rev. C76 (2007) 015801, [astro-pli/0701308]. 



-45- 



[11] D. Kharzeev and A. Zhitnitsky, Charge separation induced by P-odd bubbles in QCD matter, 
Nucl. Phys. A797 (2007) 67-79, [arXiv : 0706 . 1026]. 

[12] D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, The effects of topological charge change 
in heavy ion collisions: 'Event by event P and CP violation', Nucl. Phys. A803 (2008) 
227-253, [arXiv: 071 1.0950]. 

[13] D. T. Son and M. A. Stephanov, Axial anomaly and magnetism of nuclear and quark matter, 
Phys. Rev. D77 (2008) 014021, [arXiv : 0710 . 1084]. 

[14] J. A. Harvey, C. T. Hill, and R. J. Hill, Anomaly mediated neutrino-photon interactions at 
finite baryon density, Phys. Rev. Lett. 99 (2007) 261601, [arXiv: 0708. 1281]. 

[15] K. Fukushima. D. E. Kharzeev, and H. J. Warringa, The Chiral Magnetic Effect, Phys. Rev. 
D78 (2008) 074033, [arXiv: 0808. 3382]. 

[16] B. Link, Constraining Hadronic Superfluidity with Neutron Star Precession, Phys. Rev. Lett. 
91 (2003) 101101, [astro-ph/0302441]. 

[17] B. Link, Incompatibility of long-period neutron star precession with creeping neutron vortices, 
astro-ph/0608319. 

[18] K. B. W. Buckley, M. A. Mctlitski, and A. R. Zhitnitsky, Neutron Stars as Type I 
Superconductors, Phys. Rev. Lett. 92 (2004) 151102, [astro-pli/0308148]. 

[19] K. B. W. Buckley, M. A. Metlitski, and A. R. Zhitnitsky, Vortices and type-I 

superconductivity in neutron stars, Phys. Rev. C69 (2004) 055803, [hep-ph/0403230]. 

[20] A. Sedrakian, Type-I superconductivity and neutron star precession, Phys. Rev. D71 (2005) 
083003, [astro-ph/0408467]. 

[21] E. Babaev, A possible mechanism of effective decoupling of hadronic superfluids in precessing 
neutron stars, arXiv : 0901 . 4380. 

[22] S. Chatterjee et at, Cetting Its Kicks: A VLBA Parallax for the Hyperfast Pulsar B1508+55, 
Astrophys. J. 630 (2005) L61-L64, [astro-ph/0509031]. 

[23] J. Charbonneau, K. Hoffman, and J. Heyl, Large Pulsar Kicks from Topological Currents, 
Mon. Not. Roy. Astron. Soc. Lett. 404 (2010) L119-L123, [arXiv : 0912 . 3822]. 

[24] D. Lai, Neutron Star Kicks and Supernova Asymmetry. Cambridge University Press, 2003. 

[25] H. Spruit and E. S. Phinney, Birth kicks as the origin of pulsar rotation. Nature 393 (May, 
1998) 139-141. 

[26] C. Wang, D. Lai, and J. Han, Spin-Kick Correlation in Neutron Stars: Alignment Conditions 
and Implications, Astrophys. J. 656 (2007) 399-407, [astro-ph/0607666]. 

[27] J. M. Rankin, Further evidence for alignment of the rotation and velicity vectors in pulsars, 
Astrophys. J. 664 (2007) 443-447. 

[28] G. A. E. Wright, Pinch instabilities in magnetic stars, MNRAS 162 (1973) 339. 

[29] P. Markey and R. Tayler, The adiabatic stability of stars containing magnetic fields-ii, 
MNRAS 163 (1973) 77. 

[30] E. Flowers and M. A. Rudcrman, Evolution of pulsar magnetic fields, Astrophys. J. 215 
(1977) 302-310. 



-46- 



[31] H. C. Spruit, Origin of neutron star magnetic fields, AIP Conf. Proc. 983 (2008) 391-398, 
[arXiv:0711.3650]. 

[32] J. M. Lattimcr, M. Prakash, C. J. Pcthick, and P. Hacnscl, Direct URCA process in neutron 
stars, Phys. Rev. Lett. 66 (1991) 2701-2704. 

[33] D. Page and S. Reddy, Dense Matter in Compact Stars: Theoretical Developments and 
Observational Constraints, Ann. Rev. Nucl. Part. Sci. 56 (2006) 327-374, 
[astro-pli/0608360] . 

[34] J. N. Bahcall and R. A. Wolf, Neutron Stars. 2. Neutrino-Cooling and Observability, Phys. 
Rev. 140 (1965) B1452-B1466. 

[35] D. B. Kaplan and A. E. Nelson, Strange Goings on in Dense Nucleonic Matter, Phys. Lett. 
B175 (1986) 57-63. 

[36] N. Iwamoto, Neutrino Emmisivities and Mean Free Paths of Degenerate Quark Matter, Ann. 
of Phys. 141 (1982), no. 1 1. 

[37] O. Maxwell, G. E. Brown, D. K. Campbell, R. F. Dashen, and J. T. Manassah, Beta decay of 
pion condensates as a cooling mechanism for neutron stars, Astrophys. J. 216 (1977) 77-85. 

[38] J. N. Bahcall and R. A. Wolf, Neutron Stars. L Properties at Absolute Zero Temperature, 
Phys. Rev. 140 (1965) B1445-B1451. 

[39] G. E. Brown, K. Kubodera, D. Page, and P. Pizzochero, Strangeness Condensation and 
Cooling of Neutron Stars, Phys. Rev. D37 (1988) 2042-2046. 

[40] A. Kusenko, G. Segre. and A. Vilenkin, Neutrino transport: No asymmetry in equilibrium, 
Phys. Lett. B437 (1998) 359-361, [astro-ph/9806205]. 

[41] I. Sagert and J. SchafFner-Bielich, Pulsar kicks by anisotropic neutrino emission from quark 
matter in strong magnetic fields, arXiv: 0708. 2352. 

[42] B. L. Friman and O. V. Maxwell, Neutron Star Neutrino Emissivities, Astrophys. J. 232 
(1979) 541. 

[43] G. E. Brown, V. Thorsson, K. Kubodera, and M. Rho, A Novel mechanism for kaon 
condensation in neutron star matter, Phys. Lett. B291 (1992) 355-362. 

[44] J. E. Gunn and J. P. Ostrikcr. On the nature of pulsars. Hi. analysis of observations, 
Astrophys. J. 160 (1970) 979. 

[45] A. G. Lyne and D. R. Lorimer, High birth velocities of radio pulsars. Nature 369 (1994) 127. 

[46] B. M. S. Hansen and E. S. Phinney, The Pulsar kick velocity distribution, Mon. Not. Roy. 
Astron. Soc. 291 (1997) 569, [astro-ph/9708071]. 

[47] J. M. Cordes and D. F. Chernoff, Neutron Star Population Dynamics. H: 3D Space Velocities 
of Young Pulsars, Astrophys. J. 505 (1998) 315-338, [astro-ph/9707308]. 

[48] C. Fryer, A. Burrows, and W. Benz, Population Synthesis for Neutron Star Systems with 
Intrinsic Kicks, Astrophys. J. 496 (1998) 333, [astro-ph/9710333]. 

[49] Z. Arzoumanian, D. F. Chernoffs, and J. M. Cordes, The Velocity Distribution of Isolated 
Radio Pulsars s, Astrophys. J. 568 (2002) 289-301, [astro-ph/0106159]. 

[50] C. L. Fryer, Neutron star kicks from asymmetric collapse, Astrophys. J. 601 (2004) 
L175-L178, [astro-ph/0312265]. 



-47- 



[51] I. Sagert and J. Schaffner-Bielich, Asymmetric neutrino emission in quark matter and pulsar 
kicks, astro-pli/0612776. 

[52] I. Sagert and J. Schaffner-Bielich, Pulsar kicks by anisotropic neutrino emission from quark 
matter, J. Phys. G35 (2008) 014062, [arXiv : 0707 . 0577]. 

[53] C. Alcock, E. Farhi, and A. OUnto, Strange stars, Astrophys. J. 310 (1986) 261-272. 

[54] P. Jaikumar, S. Reddy, and A. W. Steiner, The strange star surface: A crust with nuggets, 
Phys. Rev. Lett. 96 (2006) 041101, [nucl-tli/0507055]. 

[55] M. G. Alford, K. Rajagopal, S. Reddy, and A. W. Steiner, The stability of strange star crusts 
and strangelets, Phys. Rev. D73 (2006) 114016, [liep-pli/0604134]. 

[56] S. Reddy, Novel phases at high density and their roles in the structure and evolution of 
neutron stars, Acta Phys. Polon. B33 (2002) 4101-4140, [nucl-th/0211045]. 

[57] G. G. Pavlov. O. Y. Kargaltsev, D. Sanwal, , and G. P. Garmire, Variability of the vela 
pulsar wind nebula observed with chandra. The Astrophysical Journal Letters 554 (2001), 
no. 2 L189-L192. 

[58] C. Wang, D. Lai, and J. Han. Neutron Star Kicks in Isolated and Binary Pulsars: 

Observational Constraints and Implications for Kick Mechanisms, Astrophys. J. 639 (2006) 
1007-1017, [astro-ph/0509484]. 

[59] D. Lai, D. F. ChernofF, and J. M. Cordes, Pulsar Jets: Implications for Neutron Star Kicks 
and Initial Spins, astro-ph/0007272. 

[60] M. C. Weisskopf et ai. The prospects for X-ray polarimetry and its potential use for 
understanding neutron stars, astro-ph/0611483. 

[61] K. H. Prendergast, The equilibrium of a self- gravitating incompressible fluid sphere with a 
magnetic field, i., Astrophys. J. 123 (1956) 498. 

[62] R. J. Tayler, The adiabatic stability of stars containing magnetic fields-i. toroidal fields, 
MNRAS 161 (1973) 365. 

[63] D. Page, U. Geppcrt, and M. Kueker, Cooling of Neutron Stars with Strong Toroidal 
Magnetic Fields, Astrophys. Space Sci. 308 (2007) 403-412, [astro-ph/0701442]. 

[64] M. Tinkham, Introduction to Superconductivity. McGraw-Hill Inc., 2nd ed., 1996. 

[65] L Easson and C. J. Pethick, Stress tensor of cosmic and laboratory type-ii superconductors, 
Phys. Rev. D 16 (Jul, 1977) 275-280. 

[66] A. R. Choudhuri, The Physics of Fluids and Plasmas: An Introduction for Astrophysicists. 
Cambridge University Press, 1998. 

[67] L H. Stairs, A. G. Lyne, and S. L. Shemar, Evidence for free precession in a pulsar. Nature 
406 (2000) 484-486. 

[68] T. V. Shabanova, A. G. Lyne, and J. O. Urania, Evidence for free precession in the pulsar 
bl642 - 03, The Astrophysical Journal 552 (2001), no. 1 321-325. 

[69] J. J. Matese and R. F. O'ConncU, Neutron Beta Decay in a Uniform Constant Magnetic 
Field, Phys. Rev. 180 (1969) 1289-1292. 

[70] O. F. Dorofeev, V. N. Rodionov, and L M. Ternov, Anisotropic neutrino emission in beta 
processes induced by an intense magnetic field, JETP Lett. 40 (1984) 917-920. 



-48- 



